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Abstract 

An explicit expression for the generating functional of two-flavor low-energy QCD with external 
sources in the presence of non- vanishing nucleon densities has been derived recently |fflj. Within this 
approach we derive power counting rules for the calculation of in- medium pion properties. We develop 
the so-called standard rules for residual nucleon energies of the order of the pion mass and a modified 
scheme (non-standard counting) for vanishing residual nucleon energies. We also establish the different 
scales for the range of applicability of this perturbative expansion, which are 

VQttU -0.7 GeV for the 

standard and 67r 2 /^/2mjv ~ 0.27 GeV for non-standard counting, respectively. We have performed a 
systematic analysis of n-point in-medium Green functions up to and including next-to-leading order 
when the standard rules apply. These include the in-medium contributions to quark condensates, pion 
propagators, pion masses and couplings of the axial-vector, vector and pseudoscalar currents to pions. 
In particular, we find a mass shift for negatively charged pions in heavy nuclei, AM,,.- = (18 ±5) MeV, 
that agrees with recent determinations from deeply bound pionic 207 Pb. We have also established the 
absence of in-medium renormalization in the 7r° — > 77 decay amplitude up to the same order. The 
study of 7T7T scattering requires the use of the non-standard counting and the calculation is done at 
leading order. Even at that order we establish new contributions not considered so far. We also point 
towards further possible improvements of this scheme and touch upon its relation to more conventional 
many-body approaches. 
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1 Introduction 



The pion plays a special role in nuclear and particle physics. This is related to the fact that for the light up 
and down quarks, QCD possesses an approximate chiral symmetry, i.e. it is a good first approximation 
to the theory to consider the light quarks as massless. This symmetry is not present in the ground- 
state or the particle spectrum; it is believed to be spontaneously broken down to its vectorial subgroup, 
SU(2)ixSU(2)/j — ► SU(2)i + ^, with the appearance of three (Pseudo-)Goldstone bosons which can be 
identified with the three pion states, tt^j-k . The chiral symmetry is also explicitly broken because the 
current quarks have a small mass (small compared to a typical hadronic scale of 1 GeV). The unique 
order parameter signaling this symmetry violation is the finiteness of the weak pion decay constant in the 
chiral limit, denoted /, that is / 2 7^ 0. Another order parameter often considered is the quark condensate 
in the vacuum, (0|gg|0), where |0) denotes the highly complicated vacuum. However, it is important to 
stress that, in principle, chiral symmetry could be broken even if (0|gg|0) ~ 0, as long as f 2 / 0. Due 
to Goldstone's theorem, the interactions of the pions with themselves or matter fields must vanish as 
three-momentum and energy go to zero. This in turn allows for a systematic treatment of such processes 
in the framework of an effective field theory (chiral perturbation theory, henceforth CHPT @,|U). It is 
also believed that with increasing temperature and/or density, the chiral symmetry of QCD is restored. 
While lattice studies indicate that the critical temperature is T c ~ 150 MeV, much less is known about 
the critical density. There have also been recent speculations of a very complex phase structure at high 
densities (for a recent review with many references, see It is also important to stress that lattice 
QCD applied to finite chemical potential p (density p) is only in its infancy due to the notorious sign 
problem of the Euclidean Dirac operator at p ^ (for a recent method to tackle this problem see (E[ and 
references therein). This provides another reason why it is necessary to develop an in- medium effective 
field theory. Pion properties, calculated at finite temperature and/or density, are therefore used to gain 
an understanding how these transitions are approached and in general terms to improve our knowledge 
of QCD at finite density. Another important development are the recent measurements of deeply bound 
pionic states in heavy nuclei performed at GSI HQ . These experimental data can be interpreted in terms 
of a pion mass shift due to the high nuclear density in the center of such nuclei, and these have triggered 
some calculations making use of heavy baryon chiral perturbation theory to unravel the physics behind 
this intriguing phenomenon, see e.g. refs. f|,[|. 

In this paper we will undertake a systematic study of the properties of pions and external sources 
in nuclear matter. Nuclear matter is a system of an infinite number of protons and neutrons. For 
typical densities, these interactions of the constituents of nuclear matter are strong, that is one deals 
with phenomena in the non-perturbative regime of QCD. Therefore, approximations are unavoidable, 
but one needs to be able to control these, as it is possible using effective field theory. Consequently, 
the final aim of an in-medium QCD effective field theory is to provide a systematic way of proceeding 
which allows one to estimate the errors when truncating the expansion at some order. The low energy 
effective field theory of QCD is CHPT. Indeed CHPT allows not only to tackle processes involving pions 
but as well to consider nucleons (baryons). These massive states are included as matter fields chirally 
coupled to pions and external sources. As stated, CHPT in the vacuum is a systematic way of proceeding. 
Thus there are many articles in the literature []l0|-|l3| which apply CHPT Lagrangians that are at most 
bilinear in the nucleon fields to the nuclear case in the following way: The bilinears N D N (with D 
a generic differential operator including the coupling to pions and external fields) are replaced, in the 
non-relativistic mean-field approach, by p p tvDn+ p n trD22- Here p p (p n ) are the proton (neutron) density, 
the symbol tr refers to the trace over spinor indices and the subscripts run in flavor space. Proceeding 
in this way one keeps track of the information contained in the vacuum CHPT Lagrangians, but the 
chiral counting in the medium is lost since nucleon correlations, due to the baryon propagators, are not 
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considered. In fact, such contributions can be of the same or even of lower chiral order as those terms 
kept in the mean- field approach. As shown below, they are the dominant contributions when the energy 
flowing through the baryon propagator is of the order of a nucleonic kinetic energy. However, as in the 
mean-field approach, we will not consider multi-nucleon local interactions in this paper. For attempts 



to develop effective field theories in nuclear matter without pions, see, e.g., refs. [14,TI| and references 
therein. It is also important to stress that our approach does not only encompass, but also exceeds the 
so-called low-density theorems as formulated in refs. |lq| . 

In order to go beyond the mean field approach we follow ref. || and consider the CHPT Lagrangian 
supplemented by external fields. We also use the results from ref. [|J where the in-medium contribution to 
the SU(2)xSU(2) generating functional is calculated. Consequently, the in-medium CHPT Lagrangian, 
in terms of only pions and external sources, is given. At this stage, the problem reduces to that of 
vacuum CHPT, except for the important difference that the resulting Lagrangian is non-covariant as well 
as non-local (for a general analysis of the structure of non-relativistic local effective field theories see [JT7J] ) . 
These findings are reviewed in section || and appendix |A| where the in-medium chiral Lagrangian is given 
in terms of the vacuum CHPT Lagrangian and the proton and neutron densities. The chiral expansion is 



discussed in appendix ||, while the (chiral) power counting is established in sections 2A and 2^2. We then 
apply this machinery in sections |3|— 16] to evaluate the in-medium quark condensates, pion propagation 
and masses, pion couplings to the axial-vector, vector and pseudoscalar currents and tttt scattering. In 
addition, the decay ir° — > 77 is also studied and the relevant scales of the expansion are discussed. We 
end with some conclusions in section fj]. Various technical topics are relegated to the appendices. 

2 Generating functional, effective Lagrangian and chiral counting 

For completeness, we briefly review in appendix [A] the main result of ref. [Q], repeatedly used in this 
work, where the in-medium chiral effective Lagrangian C nn is derived by integrating out the nucleon 
fields with functional methods. According to Ref. |l8| , a general chiral Lagrangian can be expanded in 
an increasing number of baryon fields, C = C n7T + Cj^ + + where ip denotes here the nucleon 

field. The results of ref. Q were obtained by truncating this series at terms bilinear in the spinor fields, 
thus neglecting the multi-nucleon Lagrangians with four or more ip fields. The reason for this is twofold. 
First, one does not know how to perform quartic Feynman path integrals exactly. Thus, the analysis 
of ref. H] cannot be extended in a straightforward way to include multi-nucleon interactions beyond a 
pure perturbative treatment of those terms already discussed in ref. [Q]. Second, and related to the latter 
point, vacuum multi-nucleon interactions are nonperturbative due to the extremely large NN scattering 
lengths and hence some kind of resummation is needed to end with an in-medium effective field theory 
of nuclear matter which also includes multi-nucleon local interactions together with pions. Nevertheless, 
the situation is encouraging due to the advances in understanding the multi-nucleon interactions in 
vacuum by the application of effective field theory to such systems. Nowadays one has two effective field 
theory schemes to tackle such problem, the original Weinberg scheme [^] (which has been made truly 
quantitative in |ll|) and the Kaplan-Savage- Weise one [p0| . We also refer to a recent and comprehensive 
review on this topic pTJ. How these advances can be applied to the nuclear medium leading to a 



satisfactory effective field theory is beyond present knowledge, see, e.g. ref. [14| for further discussions on 
this issue where the idealized case of natural interactions without pions is considered in a Fermi system. 
Nevertheless, we view it as a step forward to extend the effective field theory techniques used in vacuum 
CHPT § to the medium in the case where pions and external sources are kept, but the multi-nucleon 
local interactions are neglected. This was done in ref. jl] and here we will exploit this formalism by 
calculating systematically several in-medium Green functions within this framework. 
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In appendix |B| the expressions of the first and second order interaction operators A^ 1 ' and A^ , as 
defined in appendix |A| following ref. jy], are obtained from £ n N- The operator A is given by the difference 
between the full and the free Dirac operator (i.e. D and Dq, respectively, see eq. ( |A.3| )) and is amenable 
to a systematic chiral expansion as indicated by the superscripts "(1,2)", see eqs. QB.3| ) and ( |B.4| ), 

In sees. 2A and 2/2 we establish the chiral counting of a given diagram resulting from C^, such that 
we can determine the contributions up to a given order. We point out that one has two different counting 
schemes depending on the energy flow through the baryon propagators. 



2.1 Chiral counting: Standard case 



In this section, we establish the power counting rules of the processes that result from the in-medium 
generating functional derived in Q and discussed in appendix ffl. Our starting point is the general 
structure of an in-medium generalized non-local vertex, see fig. ffl. 






Figure 1: Diagrammatic expansion of eq. (A.l). Every thick line corresponds to the insertion of a Fermi-sea and 



each circle to the insertion of an operator T = —iA [l± — D 1 A\ . Notice that each of the three in- medium baryon 
closed loops shown in the figure is a different generalized in-medium vertex. 

As sketched in appendix |A| and detailed in jjj, each thick solid line corresponds to summing over all the 
plane wave states of the proton(neutron) Fermi-sea with three- momentum smaller than k^!\k^). In 
the following we count any Fermi momentum fcp (about ~ 2M n for nuclear saturation density po = 0.17 
fm~ 3 ) as Mjt ~ 0(p). Then each thick line, because of the three-momentum integration up to the Fermi 
momentum, is 0(p 3 ). Next, consider the non-local vacuum vertices T = —iAlI^ — D^A] -1 , see fig.§. 





Figure 2: Expansion of the non-local vacuum vertex T == —iA [I4 — D 1 A\ , where I4 is the 4x4 unit matrix. 
Every solid line corresponds to a vacuum baryon propagator iD^ 1 and each circle to the insertion of an operator 
— iA from ipDijj. 

Given the pion-nucleon Lagrangian C n N = , ipD(x)ip, the operator A is itself subject to a chiral expansion, 
A = Y ji A {i \ i > 1 (the terms in A up to 0(p 2 ) are explicitly worked out in appendix P). On the other 
hand, a soft momentum Q related to pions or external sources can be attached to each local vertex A. 
This together with the Dirac delta function of four-momentum conservation implies that the momenta 
running along the baryon propagators just differ from each other by quantities of order Q ~ 0(p). 
Since at least one Fermi-sea insertion (thick line) is required for any in-medium contribution, with on- 
shell four-momentum pp^, the four- momentum pj running along the j th free baryon propagator iDQ 1 (pj), 

#4 Note that we use the same symbol p to denote an on-shell Fermi-sea four-momentum, which is 0(1), or to indicate 
the chiral dimension of some contribution as 0(p n ) and in this case p ~ 0(p). Nevertheless, for later use, p is always 
accompanied by the symbol O or the word order and no confusion should appear. 
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is pj = p + Qj. In this way the natural chiral power counting of such propagators can be obtained by 
expanding them as, 



iD \ Pj 



. ]/+ Qj +m N . + Qj +m N 



(P + Qj) 2 -m% + ie Q) + 2Q° j E N (p) - 2Q jP + ie 
2Q v -m N + te \ 2Q v -m N ) 



Thus iD Q 1 counts as 0(p 1 ). As a result every T vertex with mp > free baryon propagators and 
(m-p + 1) A vertices (see fig. ||) scales as p VT with 

m r +l 

vp = ^2 dj — mr , (2-2) 

i=i 

where dj is the chiral dimension of the vertex Aj. Next, we will consider the chiral counting of in-medium 
non-local vertices, some of them are shown in fig. [l]. The chiral counting of an in-medium non-local vertex 
(labeled by p) that is generated through the exchange of n > 1 Fermi-sea insertions, n non-local vacuum 
r vertices, V = n + m local A vertices and m = mr^ free baryon propagators is given by 

n V V V 

d p = 3n + ^ - 4 ( n - 1) = 3n + Yj d{ ~ m ~ 4 ( n ~ ^ = Yj d{ ~ v + 4 = 4 + 12( di ~ ^ ■ ( 2 - 3 ) 

i=l i=l i=l i=l 

Here the factor 4(n — 1) originates from the four- momentum Dirac deltas attached to any T vertex. Now, 
since di > 1, any such diagram will count at least as C(p 4 ). The lowest order is obtained when only A^ 
operators with di = 1 are included. The next-to-leading order implies the inclusion of one operator 
and is 0(p 5 ).W^ 



Finally one has to take into account the exchange of pions inside or between generalized in-medium 
vertices or between them and pure pion vertices coming from the vacuum chiral Lagrangian or 
between the latter alone. Denoting by I n the total number of internal pion lines and by L x the numbers 
of pion loops, L w = I w — V p — V n + 1 where V n is the number of vertices from C wn and analogously V p 
is the number of generalized in-medium vertices. Hence a many-particle diagram originating from the 
Lag £ n N has the chiral power u, 

14 V p v-n V p V T 

v = AL 7 ,-2I n + Y,di + Y d P* = 2L n + 2 + Y,( d i " 2 ) + - 2) = 2L, + 2 + - 2) (2.4) 

i=l i=l i=l i=l i=l 

with Vt = V n + V„ and 5i is the chiral dimension of any vertex, either from C n7T or from the in-medium 
generalized vertices. We do not differentiate between di and d p i and use the symbol 5i wherever no 
confusion can arise whether a vertex comes from the free-space Lagrangian or is an in-medium generalized 



one. The counting based on eq. ( 2^4 ) is from here on referred to as the standard case. 



As a result the leading in-medium contribution begins at 0(p ), since d p > 4 from eq. (2.3) with V p = 1 



and L n = 0. In this work we will calculate the next-to-leading order in-medium corrections, 0(p ), to 

#5 In the following whenever we say that some calculation is done up to order p m , 0(p m ), we mean up-to-and-including 
contributions of order p m . It is also important to stress that the order we give always refers to the order of the terms in the 
Lagrangian (generating functional) that have been used to calculate a particular quantity. Because observables are obtained 
by differentiation of the generating functional with respect to the external sources, the leading and subsequent higher orders 
might be reduced by a common (or more) power(s) in p. This should be kept in mind throughout. 
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several re— point Green functions whenever the previous counting holds, see below. The 0(p 4 ) vacuum 
corrections to the leading non- vanishing 0(p 2 ) results can be found in ref. ||, except for the anomalous 
7T° — ► 77 decay (which only starts at fourth order in the chiral expansion since it is an anomalous process), 
and when necessary we will just give the corresponding vacuum results. As it is clear from eq. ( |2.4| ), any 
pion loop or any extra generalized in- medium vertex^ increases the chiral power by at least two more 



orders, 0(p & ). Thus the leading p— corrections, according to eqs. fl2.3| ), ( |2.4j ) are of order p 4 and p 5 and 
come from the insertion of one generalized p— vertex with only A^ 1 ' operators or with one additional 
A^ 2 \ respectively. We state once again that it is beyond the scope of the present work, and beyond the 
present knowledge, to derive a full in-medium effective field theory which includes as well nucleon-nucleon 
contact interactions. 

2.2 Chiral counting: Non-standard case 

There is one subtlety with the power counting developed so far. To understand this, we remark that the 
propagator of a nucleon in the heavy baryon formulation with four- momentum pj is + ie)- This 

corresponds to the leading non-relativistic term of the expansion given in eq. fl2.1|) which was used in 
order to determine the dimension of iD^fjjj). One caveat, already present in vacuum, arises when the 
energy is fixed and vanishes^ 7 !- In this case the heavy baryon propagator blows up although in a 
relativistic formalism used here (or in a non-relativistic one by keeping the kinetic energy) one still has 
iDq 1 (pj) = + Qj +m)/(Q| + 2Qjp + ie). The result is finite, but now iDQ 1 (pj) counts as 0(p~ 2 ) 
instead of 0(p~ 1 ) as it was used to derive eq. ( |2.4| ). To deal with this case, it is necessary to subtract 
the number of baryon lines, I B , with < kp/2niN from eq. fl2.4|) . We can write an upper bound of this 
number in terms of the number of loops, vertices and external lines. To see this note that the number 
of pion plus external source legs is greater or equal to the number of baryon lines, which can be free 
propagators as well as Fermi-sea insertions. We make use of the fact that each A vertex will involve at 
least one pion or external source, otherwise it is absorbed in the physical value of the nucleon masses, as 
explained in appendix [b|. We denote by E x the number of external legs and by I B the number of baryon 
lines minus V„, since for every generalized vertex we have at least one Fermi-sea insertion. Then we have 

E x + 2h -V p >I B . (2.5) 

Substituting in this expression I n = + Vr — 1, one obtains 

E x + 2L W + V p + 2V n - 2 > I B . (2.6) 

We can further restrict the previous inequality by subtracting 2V n on the left hand side (l.h.s.) since for 
each vacuum vertex there are at least two legs attached to it when calculating any in- medium contribution. 
Moreover, in the case of diagrams without any external source we can subtract 4V„- because any vacuum 
vertex should have at least four pion legs. Otherwise it just contributes to the vacuum renormalized pion 
propagators and associated wave function renormalizations, which can be calculated anyhow just from 
£ njT to the required accuracy. We will denote by a) the first case and by b) the more specific second one. 
Then we can write: 

a) E x + 2Ar + V p - 2 > I B , 

b) E x + 2L n + V p -2- 2V n > I B ■ (2.7) 



#6 Henceforth we will independently use the words in-medium or the symbol p— to denote any contribution due to the 
finite nucleon density of the ground state. The symbol p— should not be confused with the p meson, which we always denote 
as p(770). 



#7 The expansion of eq. (2.1) is of quantities Q 2 ~ M% ~ kp over 2Q°mjv and then, when is of the order k F /2mN, it 
breaks down. 
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We now indicate by n m the number of A vertices with m attached legs. Then for each of such vertices 



m — 1 lines are non-baryonic and can be removed from the l.h.s of eqs. (2/7). In the same way, we call p r< 
the number of vacuum vertices from C nn with m legs. For the case a) we have already subtracted — 2V 7T 
and for the case b) —4V n . Thus we can further remove m — 2 for a) and m — 4 for b). Putting all this 
together, we arrive at the equalities: 



a) E x + 2L n + V p - 2 - ) (m- l)n m -}(m- 2)p m = I B , 

m>2 m>3 

b) E x + 2L w + Vp-2- 2V n - (rn - l)n m - £ (m - 4)p m = I B . (21 



m>2 m>5 



Notice that the second of these equations is just a rewriting of the first one for the special case with only 
pion legs. In summary, after the elimination of L n and the subtraction of I B , the counting index v given 
in eq. (|2.4| ) changes to 



a) v = 4 + J2( S i " 3 ) + - 2) - ^ + X) ( m - + X ( m ~ 2 ) Pm + ( /s " 7 ^ ' 

i=l i=l m>2 m>3 

b) i/ = 4 + ^(Si - 3) + S i ~ E * + X) ( m - l)"m + X ( m - + ( /b - J b) • ( 2 - 9 ) 

i=l i=l m>2 m>5 

When 14V = 0, we adopt the convention that there is no contribution from the second and fourth sums on 
the right hand side (r.h.s) of the last equation. The interesting point of these equations is that for a given 
number of external legs E x the chiral dimension v is bounded since b% > 4 for a generalized in-medium 
vertex as shown above and b% > 2 for any vertex from Cmr and by definition Ib> !%■ Note also that the 
inclusion of an extra generalized p— vertex increases v at least by 1. Those loops that arise because of 
pion lines that are exchanged inside only one generalized p— vertex will increase Is — 1% by the number 



of involved baryon propagators in the loops because of the absence of the so called pinch singularity [18|. 
However, when two or more p— vertices are involved in the pion loops this is no longer the case. We will 
apply these counting rules when discussing in-medium pion condensation and 4ir scattering. In the latter 
case, we will just calculate the leading contribution 0(p 3 ). The counting based on eq. ( |2.9D is from here 
on referred to as the non-standard case. 

In the sections |3]-||, we will use the simpler counting rules, see eq. ([2.4|), to determine the diagrams 
to be considered up to 0(p 5 ), with the exception of those cases cases where a small energy denominator 
occurs, <3° ~ kp/2mN, and where we will turn to eq. (|2.9| ). This will always be indicated in the text at 
the appropriate places. 

Finally, let us stress again the difference in the power counting for pion properties in the medium and 
in the vacuum. While in vacuum the leading and next-to- leading order is 0(p 2 ) and 0(p 4 ), respectively, 
in the medium and for the standard case, these contributions start at 0(p 4 ) and 0(p 5 ), respectively. 

3 In-medium quark condensates 



From the generating functional Z(v, a, s,p) given in eq. ( |A.l ), the quark condensate is obtained by partial 
functional differentiation, 

(n\qiqj\n) = -- — —Z(v,a,s,p)\ v ^ p= o , (3.1) 

OSij [X) 
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with qi = u and q 2 = d quarks. Up to 0(p 5 ), because of the absence of pion loops, the in-medium 
contributions can be calculated directly by equating the generating functional Z with the in-medium 



action j dxC 7Tn (x) and the vertex can be identified by doing the partial derivative 



dsij(x) 



X). 



The 



resulting p— diagram is depicted in fig. ^ with one insertion from since the scalar source s(x) is not 
present in (see eqs. QB,3j ) and ( p3.4j ) ) . 




Figure 3: Unique diagram that contributes to the calculation of the in-medium quark condensates up to 0(p 5 ). 
As usual the thick line represents a Fermi-sea insertion, the circle is an vertex and the wavy line an external 
scalar source s(x). 



We denote by p p and p n the proton and neutron densities, in order. These are given by 



Pp 



f dp 



(2vr) J 



3vr 2 



Pn 



k f dp (kp) 



(2vr) J 



3tt 2 



(3.2) 



and we also introduce the quantities p = (p p + p n )/2, the isospin symmetric nuclear density, and p = 
(p p — p n )/2, the isospin asymmetric nuclear density. It is easy to see from eqs. ( A.l ) and ( p.4| ) that the 
in-medium contributions to the action j dx C which contain the scalar sources (the ones relevant for the 
calculation of the quark condensates) are: 



dx 



dp 



(2n) 3 2E( y p) 



Tr 



{ Cl 4Bo(s(x)) + c 5 [4B s(x) - 2B (s(x))}} n(p)(]/+m N ) 



dx 



4B p(2 Cl - c 5 )(s n (x) + s 22 (x)) + 4B c 5 ( Sll (x)p p + s 22 (x)p n ) + 0(p 6 ) , (3.3) 



where the usual SU(2) matrix notation is applied for the scalar source. By differentiating with respect 
to sn(x) and s 22 (x) one obtains, respectively, the up- and the down-quark condensates at 0(p 5 ): 



{n\uu\fl) 

(n\dd\n) 



(uu), 
{dd) y 



1 



2a A 4c 5 _' 
pMf + p P 
2a A 4c 5 _" 



(3.4) 



with the pion-nucleon sigma-term a = —Ac\M% at 0(p 2 ) P^,p3[. Furthermore, the subscript vac refers 
to the vacuum value of the corresponding quantity up to 0(p 4 ) as given in ref. Q, which at lowest order 
is simply -B p. 

Keeping terms up to the same accuracy,^ the isoscalar sum (Q\uu+dd\ft) = (uu+dd) vac [l — 2a pj f 2 M%\ 
only depends on the total nuclear density 2p = p p + p n and agrees with the expression, valid up to linear 

#8 In the following, when the previous sentence is not explicitly stated, the context will clarify when an equality is exact 
or accurate up to higher order corrections. 
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order in the density, given in ref. |M| under the use of the Hellmann-Feynman theorem (we are using a 
chiral perturbative approach and the agreement is only valid strictly when calculating the cr-term up to 
0(p 2 ) as noted above). For a more intuitive derivation, see ref. p5| . This coincidence is not surprising 
since, in the non-relativistic limit, the diagram of fig. |3| just counts the number of nucleon states times 
the corresponding vertex for protons and neutrons separately. If one assumes an expansion in density 
on intuitive grounds this has to be the leading result. This is indeed the argument used in ref. [25|. 
Note that deviations from the relativistic limit in fig.||] increase the counting by two orders. As a new 
aspect with respect to previous works 12, 24], we also consider strong isospin breaking represented by 
the term with C5 in eq. (|3.4j ). Other works that also consider the calculation of the quark condensates 
in symmetric nuclear matter are refs. [11, 12| within a mean-field approach and refs. [26 29| using the 
Nambu-Jona-Lasinio model. Note also that the quark condensates are essential inputs of in-medium 
QCD sum rules techniques [^, 30 1 . 

For the case of symmetric nuclear matter, all the previous investigations led to dropping quark conden- 
sates with increasing density. Specifically, using the value for the cr-term, a = 44 ± 8 ± 7 MeV one 
obtains a relative reduction of a (35 ±9)% with respect to its vacuum value for nuclear matter density. 
This reduction of the quark condensate with density is traditionally considered as a clear indication of 
chiral symmetry restoration in the nuclear medium p6| . At this point it is worthwhile to stress that 
the quark condensate is not a unique order parameter of chiral symmetry breaking. Even if it vanishes, 
chiral symmetry is still spontaneously broken as long as f 2 7^ 0. We will come back to this point below 
when evaluating the in-medium weak decay couplings of the pion. We also note from eq. ([TJ) that the 
correction becomes unity when p ~ 2.8po corresponding to a Fermi-momentum kp ~ 400 MeV. However, 
this value should also be considered indicative since then the correction is equal to the leading term. 

Consider now the isospin breaking contribution. From ref. [23] we take the value C5 = —0.09 ± 0.01 
GeV^ 1 . For realistic cases p n > p p and then since C5 < 0, the in-medium dd condensate is smaller than 



the uu one. Using (uu) 



(dd), 



in the limit m u = ||, the pertinent difference is: 



4c 5 



Pn - Pp 



(n\uu - dd\n) = (qq) vac -prP ^ (qq)vac 0.028 

/ Po 

which is very suppressed as compared to the density dependence of the sum given by 



(3.5) 



(tt\uu + dd\tt) ~ (qq)^ 



1 - 0.35 



Pn + Pp 
Po 



(3.6) 



because — 2cs, see. eq. (|3.4|), is more than one order of magnitude smaller than cr/M 2 ~ 2.3 GeV . Hence 



the presence of the terms in eq. (3.4) proportional to C5 does not alter appreciably the reduction of the 
quark condensates at finite density present in the symmetric nuclear matter case studied in refs. [12,24, 2^]. 



4 Two-point Green functions 



We now consider the propagation of pions in the medium. We will first derive their equations of motion 
from which one can read off the in-medium pion mass. We will also discuss pion-condensation illustrating 
the use of eq. ( |2.9[ ). In addition the new scales for the in-medium chiral expansion are also pointed out. 
The p— wave function renormalization of pions needed to satisfy the quantum canonical commutation 
relations is derived in appendix Finally, we will also study the couplings of pions to the axial-vector 
and pseudoscalar currents in the medium and their relation as dictated by a QCD Ward identity. 



9 



4.1 Equations of motion 

In fig. |I| we show those diagrams that give the in-medium corrections of the pion propagator up to 0(p 5 ). 
As pointed out in ref. ||], the p— vertices have analogous properties to those of the standard local vertices 
(i.e. the field theoretical vertices from C n7T ) in the computation of the numerical factors accompanying 
the exchange of pion lines in a given diagram. With respect to fig. || this manifests itself in the fact that 
together with the momentum Q flowing from the left to the right, as in particular shown in fig. |I]b, one 
has also to consider the inverse process, Q — > —Q. This is immediately accomplished when working in 
configuration space. 



a) b) c) d) 

Figure 4: Diagrams that contribute to the in-medium propagation of pions up to 0(p 5 ). Figs. a), b) and c) are 
0(p 4 ) and fig.d) is 0(p 5 ). The crossed diagram of b), obtained by replacing Q — ► — Q, is not shown. The energy 
through the pion lines (dashed) is ~ M T . 

Denoting by Q the pion four-momentum and by p the one running in the (lower) Fermi-sea insertion 
(p 2 = m 2 N ) in each diagram of fig. ^j, one can easily see that the imaginary parts of the diagrams in this 
figure vanish for Q° ~ M n 3> Q 2 /(2mjv)- Note that the presence of the pole of the baryon propagating 
with four-momentum p — Q in fig. || implies the vanishing of the denominator (p — Q) 2 — fn% = Q 2 — 2pQ. 
This occurs only when Q° ~ Q 2 /2m^ which is not the case for Q 2 ~ M 2 , i.e. for the type of pion 
momenta considered in most parts of this work. The same reasoning can be applied to fig. Denoting 
by q the four-momentum of the second Fermi-sea insertion we then have (p — Q) 2 = q 2 = m 2 N which is 
once again the same equation as the one required to pick the pole in fig.[|b. In fact, the imaginary parts 
of fig. ||b and fig. ||c are closely linked and tend to cancel each other. Remembering the Feynman rules 
of appendix |Aj a second Fermi-sea insertion implies an extra factor of 

1 f kp da 

' 1 -W+m N ), (4.1) 



2 J (27r) 3 2£(<7) 

as compared to the case with only one Fermi-sea insertion. On the other hand, putting a baryon on-shell 
in fig. [|b implies 

dq </ -\~TYl n If c?q 



I 



(2vr) 4 q 2 _ ^ + %e 2 J (27r) 3 2E[q) V ^ +m ^ ~ ■ 

with the four-momentum q = (E(q), — q). The second term on the r.h.s. of this equation comes from the 
anti-baryon pole and does not contribute because of the presence of an extra delta function 6(q° —p° — Q®) 
which cannot be satisfied since in this case q° = —E(q) and p° is an on-shell baryon momentum coming 
from the sum over the states in the Fermi-sea as already discussed. Then only the first term on the r.h.s. 
of eq. (^[^) is possible and by splitting the three-momentum integral as f°° <iq = J kp (iq + dq, we 
see that the integration up to &f cancels against eq. ( [4.1|) , leaving only dq. This implies that the 
absorptive parts, arising by exciting one baryon on-shell, begin to contribute only when the baryon-three- 
momentum is larger than &f (particle excitation), as one could naturally expect since all the other states 
in the Fermi-sea are already occupied. 
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Introducing the density matrix p = diag(/5 p , p n ), the in-medium contributions to the quadratic pionic 
term of C, SC^, can be obtained from eq. ( |A.1[ ) (as just discussed, the last term of eq. ( A.l ) with two 
Fermi-sea insertions does not contribute here). We then obtain 

i J dx5C H {x) = J dx{p[4>,^\) + J dx^pi - 2c l B ( ( f> 2 M) + ^02(H>) + -cg^W) 

- ^c 5 B o {4> 2 M+M0 2 + 2</>M<fi) +c 5j B o (0 2 7W)|) (4.3) 

dxdy e^ x -^Tr{ 7 ^ 5 d^(x)DQ\x,y)Yl5duHyW +m N )n(p)\ , 



4/ 2 J (2tt) 3 2E( P ) 

where the Cartesian coordinates for the pion fields are defined by eft = Ya=i with 4> given in eq. (|B.2| ) 

o 

and Ti are the usual Pauli matrices. In the last expression we have set 9a to its physical value qa = 1.267 
since corrections will be of higher orders than the next-to- leading one, qa =9 a {1 + 0(rh)} p2| , and for 
the same reason we take in the following m p = m n = mjy = (m n + m p )/2 ~ 939 MeV as mentioned in 
appendix |A]. 

The in-medium contribution to the equations of motion are obtained by performing the differentiation 

dQ 

9aP (Q 2 ) 2 m^ -9aP{„ Q 2 \ , -Pr> ^ 16£? ° 



5 f dyC^j '54>k{x). Working in momentum-space, 4>k{x) = f , 2 Ri e l ® x 4>k{Q)i one obtains: 



4>{Q)k - i^jf (Qo ~ ^\ e k j3(p(Q)j + ij2Qo£kj3(f>(Q)j ~ ^^pmc!(j)(Q) k 



2f 2 m N Ql 

f 2 (^Qo + ^Q 2 )HQ)k-j 



+ {c 2 Q 2 o + c 3 Q 2 ) <KQ) k - pmc 5 5 3k cf>(Q) 3 . (4.4) 



Diagonalizing these equations by working in terms of the charged fields 7r + = (<p\ — i(\>2)j\[2 and ir~ = 
(01 + i<p2) /\/2; together with 7r° = 03 and adding the vacuum piece, we finally have the following spectral 
relations between the energy Qq = oj and the three-momentum Q for on-shell in-medium pions: 



7T 



-M 2 -£-c 5 = , 
j z m 



/ rn N f z 



9\P (Q 2 ^ 2 

2f 2 m,N oj 2 



7T 



. 2 -M 2 + (l + Cl ^ + ^ 2 U + <*-^) 



PJ P 



+ P^ P 



Q 2 1 + ^C 3 



4 P Cj 9 2 aP 
P rn N p 



9\P (Q 2 ) 2 
2prriN oj 2 



7T 



^-m 2 (i +c1 8 4)+% 2 (c 2+C3 -A-) 



pj p 



A 

Stun J 



pu 4 + p U 



Q 2 l + ^C 3 



4 P Cj 9aP 
P rn N p 



9\P (Q 2 ) 2 

2/ 2 mjv 

(4.5) 
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with M 2 = 2mBo the lowest order CHPT pion masses and M n o, M n + the vacuum pion masses at 
C(p 4 ) f3|. We will take all of them to coincide with the corresponding physical masses since differences 
will be 0(p 6 ). 



4.2 Pion dispersion and condensation 



The previous expressions largely simplify in the case of symmetric nuclear matter p = in the chiral 
limit, 



2 2 L 4 P 
w =Q { p° 2 

The pion velocity in-medium for this case is given by: 



duj 



i 



2p 

P 



72 C 2' 



(4.6) 



(4.7) 



Since it must be smaller than the velocity of light 17, 33|, this imposes the constraint ci > which is 
satisfied by the actual value of this constant, C2 = 3.2±0.25 GeV -1 |p4[| . Imposing v > 0, then C2 < f 2 /2p. 
Taking the previous value for C2 gives p < 2pQ. This clearly indicates that in-medium CHPT, just by 
considering irN interactions with the actual values of the c\ counterterms, can only be applied to densities 
< po, as otherwise the corrections will be larger than 50%. For po (nuclear matter saturation density), 
kp — 270MeV~ 2Ai 7r . Even for vacuum ttN scattering, such a value for the running three- momentum 



is on edge of applicability of the theory [34,35]. However, one also has to take into account that due to 
the large variation in the values of the c« and combinations of them, important differences with respect 



to the actual scale for a specific process are expected; see section O for a discussion about the scales at 
which the chiral expansion breaks down. 

We are now interested in finding solutions to eqs. Q4.5|) with vanishing energy so that the ground state 



becomes unstable against the spontaneous creation of pion modes. A close look to eqs. (4J3) indicates the 
presence of possible solutions in the antisymmetric matter case, with Q, 2 /u> = constant for u, Q — ► 0. To 
be more specific, let us consider eq. (4.5) for negative pions ir~ . Then for u, Q — > and Q 2 fu 
has the equation: 

8C1 ^> <■ (4.8) 



D one 



9aP 
2pm N 



D 2 + 9 ^-D + Ml + 



P 



1 + 



P 







with the solution: 



D 



Pm N 
9aP 



9aP 




Pm N 



M 2 



1 + 



8cip 



(4.9) 



In order to have an oscillatory mode and not a diffuse one, the following inequality has to be fulfilled, 



P> 



4/ 2 



g A m N \p 



M 2 



1 + 



(4.10) 



When the previous equality holds, D = —m^p/p. This implies that for the tt + , where D = m^p/p, 
the energy is negative so that it is energetically favorable to excite tt + modes, i.e. we have a n + boson 
condensate. Because of this condensate, the U(l) flavor symmetry, which is still present in eq. (|A.l| ), is 
spontaneously broken, leading to one Goldstone boson and two heavy modes because we started with 
three independent pion modes. See ref. for a similar discussion in the case of QCD at finite isospin 
density pi in the idealized case of pb = 0. Note also that for the neutral pion there is no such solution. 
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For \p\ -C p, it is justified to use eqs. ([4.5D because the expansion of the baryon propagators in eq. ([O) 
still holds since then D <C Im^. However, when \p\ ~ p the remarks discussed at the end of sec. |2.l| 
apply and a more careful treatment is necessary. Still, the prior situation is approximately realized in 
heavy nuclei, where for N/Z = 1.5 one has p/p = —1/5. Thus we can apply eq. (4.10) with the result: 



so that: 



P 



> 8 P0 . (4.12) 



The previous inequality implies: 

p > (1.6 ± 0.2)p , (4.13) 

where we have used c\ = — (0.81 ± 0.12) GeV -1 |pS|| . The resulting lower limit for the density in eq. (|4.13j ) 
is larger than the total density in heavy nuclei, p ~ pq. Therefore, these oscillatory pion modes in heavy 
nuclei do not occur. 

For neutron star matter, p = —p, there is a reduction by a factor of 25 on the l.h.s. of eq. ( 4.1 1| ) so 
that: 



'-^ M -H 1 + 7^°M 1 + 7^' (4 ' 14> 



which amounts to a drastic change. In fact, for this case, eq. (2.1) cannot be applied and one has to take 
care of the modification of the counting by making use of eq. (2.9) case b) since here we have only pion 
lines. 

In appendix |C] we discuss in detail the use of the non-standard counting rules applied to the problem 
of the pion propagation in the medium. These results will be used in the next section when determining 
the appropriate in-medium scale in our framework. 

4.3 In-medium breakdown scales 

In this section, we will deduce the in-medium scales at which the chiral expansion breaks down. To 
do so let us consider eq. ( A.l| ) for the generating functional or better its more general, non-perturbative 
expression given in ref. [Jl|. There, the in-medium contribution to the chiral Lagrangian is given by: 

/ dxdy / feF (27)feb) e " P(X_y)Tr{(]/+mAr)l0g - F} ' (415) 



with log J 7 being a pure interaction operator whose perturbative expansion leads to eq. ( |A.1| ) (for defini- 
tions, see ref. |TJ). In the limit kp — > 0, the leading non-relativistic contribution, if not vanishing, comes 
from ]/ +uin — > 2 itin and the above expression can be written as 

J dxdyTr{2m N log^} J** —J*. = J d ^ (J dyTr{log^}) , (4.16) 

involving the numerical factor 1 / 6n 2 and establishing that the in- medium corrections must start at order 
k F oc p. Taking into account the perturbative expansion of log given in eq. ( [All) and noting that 
the local operators A have dimension of mass, see e.g. eqs. ( |B.3| ) and flB.4j ), one concludes that any term 
arising from eq. (|4.16| ) with m pion fields and n momenta, n > 1, counts as: 



(4.17) 



13 



with Q n involving both derivatives and quark mass insertions and A refers to a typical hadronic mass 
~ 1 GeV. On the other hand from C nn one typically has, n>2: 



f'm-2 \n-2 ' 



To compare eqs. ( 4.17 ) and ( 4,1§| ) it is convenient to rewrite the former as: 



kp 



6vr 2 / 2 f m ~ 2 A n - 



(4.18) 



(4.19) 



Since in our counting kp ~ Q we should consider n —* n — 1, n > 2, in eq. ( 4.19 ) so that the local in- 
medium contributions are suppressed by k F /{Qirf 2 ) as compared to the vacuum ones. It is also worthwhile 
to stress that while the leading in- medium contribution scales as kp ~ p when kp — ► 0, eq. ( [4.16 ), the 
suppression with respect to a large scale surviving in the chiral limit, much larger than pion mass or 
momenta, involves the second and not the third order in &f, i.e. the suppression is non-analytic in the 
nuclear density. Note that this scale, v6/ — 700 MeV, can also be inferred from the S-wave Weinberg- 
Tomozawa contribution (the last term in eqs.( [4.5D for the charged pions). Under P-wave interaction, see 
the last but one term in eqs. (4J5) for the charged pions, the scale is reduced to y/Qf jgA — 570 MeV, 
a little bit larger than twice the Fermi momentum for nuclear matter saturation density. It is roughly 
a factor of two smaller than the commonly used value 47r_/' 7r typically quoted in vacuum CHPT for the 
pion sector. This is similar to the flavor dependence ~ 1/Nf pointed out in the context of extended 
technicolor approaches to electroweak symmetry breaking |^6| . 

The main problem when Q® ~ Q 2 /2mj^, i.e. the non-standard scenario for the counting rules, comes 
from the further reduction in the scale A of the nuclear perturbative expansion pertinent to such cases. 



In fact, instead of having the factor 2rri]yQ as in eq. (2.1), one needs to consider the full expression since 
Q° counts now as a quantity of order k F /2mjy in the denominator, kp ~ 0(p). As a consequence, if 
something counted before as k F /A 2 , it now scales as 2m^kpl A 2 . This implies that A —* A 2 /2mjv — 170 
MeV instead of A ~ 570 MeV, that is, a reduction by roughly a factor of three. For instance, considering 
eq. ( C.l ) for the propagation of pions in the case of u = we have from fig. [|a, which in fact corresponds 
to figs.0b,c: 



2m N g 2 A Q A 
f 2 



kp 



+ 



4k 2 p 



log 



2k F + |Q| 



32vr 2 |Q| b |2/c F -|Q| 



Taking now |Q| <C kp, the previous expression reduces to: 



Q 



2 m N g 2 A k F 
^ 2 f 2 



(4.20) 



(4.21) 



with the scale A' = Air 2 f%/2m N g 2 A ~ 120 MeV, of the size we have estimated just above. As a result, 
while for the standard counting scenario one has a chiral expansion to be applied when kp < 2M n , 
corresponding to nuclear saturation density, when Q® ~ Q 2 j2m^ the perturbative expansion can only be 

applied up to densities around ~ (170 MeV/A;^) 3 ~ /?o/4, where k F ~ 268 MeV is the Fermi momentum 
at nuclear matter saturation density. This modification of the scale is simply lost in the mean-field 
approaches to nuclear matter. However, as we have seen in sec. 4.2 and appendix^], this scale is crucial 
for the study of pion condensation. 

For the case of the K~ condensation in nuclear matter studied in [1C] under the use of the mean-field 
approach, the absence of baryons with strangeness +1 implies that an analogous diagram to that of fig. ^a 
is not involved. 
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a) 




6 



c} 



Ya |,| )-c-{a" i Wa i 




d) 



e) 




Figure 5: Diagrams that contribute to the in- medium propagation of pions up to 0(p 4 ) when Q° ~ kp/2m,N, 
eq. ( |2.9| ) case b). Diagram a) is 0(p 3 ) and the rest are 0{p A ). Other diagrams obtained by exchanging free baryon 
propagators by Fermi-sea insertions or by exchanging lines and or the position of loops are not shown. 



The previous problem is similar to that of studying nucleon-nucleon interactions in the vacuum where 
one also has this situation [18] due to the appearance of the large 2m n ~ 0(1) factor instead of an 
0(p) one. The solution advocated in this case is to perform a resummation of diagrams by solving the 
Lippmann-Schwinger equation [O] and resumming those diagrams with two nucleons intermediate states. 
It is clear that one should proceed in a similar way here by resumming an infinite subset of diagrams, 
including as well nucleon contact interactions, as, on the other hand, is known for a long time in many- 
body problems [43,44]. More specifically, including extra p— vertices repeatedly iterated as in fig.|5^, with 
any of them enhanced by the appearance of the large factor 2mj\i/Q, amounts to the RPA approximation 
for the in-medium pion propagator. 



4.4 In-medium pion mass 

There has been recent interest in calculating the pion masses in the medium using heavy baryon chiral 
perturbation theory (HBCHPT) Both works include one pion loop graphs, C(p 6 ), and while ref. || 

is restricted to the nuclear symmetric case, ref. |9j also considered non-symmetric nuclear matter although, 
as shown below, not all the effects of order 0(p 6 ) are included. However, the subject has already a long 
history and similar diagrams than those calculated in refs. MM were already considered in earlier (chiral) 



models [37], and the pole position of the pion propagator in nuclear matter has been discussed in the 



context of pion or kaon condensation, see e.g. [38|. The calculation of the in-medium pion mass implies 



to impose in eqs. (|4.5|) the condition Q = and then to solve for to = M n , where the tilde denotes 
in-medium pion mass (afterwards, when discussing a QCD Ward identity and pion condensation, the full 
spectral relations will turn out to be essential). From eqs. fl4.5|), the equations to determine M n are 



tt : 



7T 



MJ, - Ml, (, + Cl f) + (* + c 3 - 4) - Mfp s - , 

M> + -M^ 1+ ^) + ^ C2 + C3 -^)-^0, 
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7T 



M 



M 



1 + Cl 



8/5 \ 4p 



C2 + C 3 



A. 

8mN J 



\ + PM JL 



(4.22) 



Up to the order we are considering here it is correct to substitute M n — * M n in all the expressions above 
except when M% appears alone in the beginning of these equations. We will specialize our study to the 
case of Pb, a heavy nucleus, where deeply bound ir~ states have been detected |],[7j with a shift in the 
effective in-medium tt~ mass AM n - = 23 — 27 MeV [f?J under the use of pion-nucleus optical potentials. 
To facilitate the comparison with ref. || we take the same values of k p = 246.7 MeV and k n = 282.4 MeV 
corresponding to almost nuclear matter density saturation and N/Z = 1.5. These are typical values for 
heavy nuclei. 

As values for the next-to-leading order pion-nucleon counterterms we take from ref. [39| c\ = —0.81 ± 
0.12 GeV" 1 , C3 = —4.70 ±1.16 GeV" 1 , whereas C2 is essentially undetermined by this analysis. The 
value of C3 has been also determined in ref. |^0| in good agreement with the previous value. Finally from 
ref. [34, 35 1 we take C2 = 3.2 ± 0.25 GeV -1 . A look at eqs. ( |4.22| ) reveals that the previous formulae only 
depend on the combination C2 + C3 which can be pinned down more precisely from the experimental value 

Of T + (M n ): 



0.045 ± 0.088 fm gy. Following ref. H one has: 



T+(M n ) = -f 



-4ci + 2c 2 



+ 2c 3 + 0(p 3 



(4.23) 



Considering the previous value for c\ one finally has C2 + C3 = —1.46 ±0.26 GeV 1 with a central value in 



agreement with the values for C2 and C3 quoted above. In fact, taking the value of C2 from [34] it follows 
that C3 = — 4.66±0.36 GeV -1 , which is the value we will finally take. It is more precise than that of ref. [39 
but with the same central value. From the previous equation we can also determine the combination 
— 2ci + C2 + C3 = 0.16 ± 0.10 GeV -1 . This result is consistent within errors with the one obtained from 
analyzing pion-deuteron scattering in chiral perturbation theory, — 2ci±C2±C3 = — 0.04±0.02 GeV -1 [42|. 
The difference in the central value reflects an 0(p 3 ) contribution to T + (M 7r ) ||23[| . Nevertheless, when 
this combination of counterterms appears we will directly take the empirical value of T + (M 7T ). 





Exact (AM) 


Exact (M) 


Perturbative 


AM ff - 


18 ± 5 MeV 


18 ± 16 MeV 


8.2 ± 2 MeV 


AM 7r + 


-12 ± 4 MeV 


-12 ± 13 MeV 


-6.5 ± 2 MeV 


AM n o 


2 ± 4 MeV 


2 ± 7.0 MeV 


1.1 ± 2 MeV 



Table 1: Shift between the in- medium pion effective masses and the vacuum 
ones. The label Exact (AM) refers to solve eqs. ( 4.22j ) exactly in terms of 



5uj' 2 = Ml — Ml and analogously for Exact (M) but directly for M v . Finally, 
Perturbative means to solve the eqs. ( [4.22| ) by replacing M n by M n in all the 
terms with density dependence. 

In the third column of table [l] we show our numerical results for 207 Pb in the perturbative case when 
substituting in eq. ( [4.22 ) M n — * M n in all the terms with density dependence. As we can see the resulting 



shift in the mass of the tt~ is much smaller than the experimental one, AM n - = 23 — 27 MeV. However, 
looking closely at eq. ( |4.22| ) one realizes the presence of large in-medium corrections, of the order of 
50%. These originate from the terms within the round brackets, namely 1 + 4p(c2 + C3 — g\/8ni]y)/ f 2 — 
1 ± 8pci/ f 2 ~ 0.5 for nuclear matter saturation density. Thus doubts arise with respect to the accuracy 
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of eqs. ( 4.22| ) indicating that higher order corrections will give non- negligible contributions. However, 
this also indicates that if one gives credit to eqs. ( 4.22 ) for such densities, one has to solve them exactly. 
Otherwise one is ignoring large corrections to the solutions of eqs. ( |4.22| ) when solving them in a pure 
perturbative way. Note that if 5 = —0.5 such that 1 + 5 = 0.5, then the deviation from one for ratios 
like 1/(1 + <5), appearing when solving eqs. (422), is a 100% effect. On the other hand, solving eqs. ( 4.22j ) 
in terms of M n , as they are written, one ends with the results collected in the second column of table [l]. 
When comparing columns two and three two important points are easily noted: 1) The sizeable difference, 
more than a factor of two, in the central value of the shift of the 7r~ mass when comparing the full and 
perturbative solutions and 2) the errors for the second column are so large that they prevent one from 
obtaining any definitive conclusion about the mass shift. This situation can be improved if one rewrites 
eqs. ( [4.22 ) in terms of 5uj 2 defined by the change of variables M 2 = M 2 + 5uj 2 and then solving for 8iv 2 
in the ir~ , ir + and 7r° cases. In this way one is sensitive to T + {M 1T )p and to C2 + C3. Furthermore, 
the square of the vacuum pion mass is removed in the equations, which sizably decreases the sensitivity 
with respect to C2 + C3. Proceeding in this way one obtains the numbers in the first column of table |l| 
that constitute our final results. Within errors, these are compatible with the experimental interval of 
AM n - = 23 — 27MeV. 

Higher order corrections to this result have to be evaluated but not using M n = M w as it was done 
in ref. ||. Indeed, we also indicate here that the corrections to M 2 in an asymmetric nuclear medium, 
because of the Weinberg term (the last terms in eqs. ( 4.22j ) for the charged pions), are C(p 4 ). Hence if one 
is interested in an C(p 6 ) calculation as in ref. Q, one should take the value at C(p 4 ) for M v and use this 
one to determine once again the Weinberg term. If this is not done, one loses 0(p 6 ) contributions which 
in fact were not taken into account in ref. Q because there the in-medium self-energy was calculated 
with Q = (M n , 0). However, such kind of corrections are numerically small as they contribute just ~ 0.5 
MeV to AM,-. 

To conclude, our results in the first column in table [l] tentatively indicate that in-medium CHPT 
accounts for most of the required shift in the measured ir~ mass at finite density from recent experiments 
on deeply bound pionic atoms 0. The main contribution to this shift, around 16 MeV, results from the 
division of the Weinberg term by the factor 1 + 4/3(c2 + C3 — g\/8mj^)/ f 2 . The latter factor corresponds 
to the wave-function renormalization of the in-medium pions in symmetric nuclear matter at threshold, 
c.f. eq. ( p. 15 ), which for nuclear saturation density is ~ 0.5. This last point, related to the propagation 
of pions in the medium, was not accounted for in ref. ||. Notice as well that the first term in brackets in 
eqs. ( [4.22 ) is also around 0.5, such that, when divided by the wave-function renormalization constant, it 
is slightly greater than one and it further increases M n - by the small amount of 2 MeV, the same value 
appearing in the M% case, giving finally the 18 MeV of increase for the tt~ case reported in table |l|. 



4.5 Axial-vector and pseudo-scalar pion couplings 

The in-medium axial- vector current A l ^{x) = 9(^)7^75 (t* /2)q(x) is given at the tree-level by a calcu- 
lation of the functional derivative of the classical action j dyC[a,v, s,p] with respect to a^(x) where 
a{x) = Y2i a li.( x ) Tl /2 an d where, in the end, the limit of vanishing external sources is taken. In a dia- 
grammatic language and for the in-medium contributions, this amounts to consider analogous diagrams 
to those of fig.|| but with one external pion line substituted by an axial- vector source. In addition, one 
needs to consider the pion wave function renormalization. This deserves special attention in our case be- 
cause the non-relativistic and non-local character of the in- medium chiral Lagrangian deduced in ref. Q , 
as also discussed in appendix [A|, eq. ( |A.1| ). Thus standard prescriptions to obtain the wave function 
renormalization are not appropriate here and a detailed calculation of the latter is given in appendix [D[ 
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We also introduce the axial-vector current J^ 2 



4 1 



d"lfi75U and its hermitian conjugate 
(Jj^jt. Furthermore, because of the breaking of covariance due to the presence of the nuclear medium, it 
is convenient to separate between temporal and spatial couplings of the pions to the axial- vector currents. 
In this way we have: 
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,(4.24) 



where f n = 92.4 MeV is the vacuum weak pion decay constant. Notice that the terms in curly brackets for 
the temporal components correspond to the square root of the wave function renormalization, -Z(Q) 1 / 2 , 
eq. (D.14), only in the case of symmetric nuclear matter. Furthermore, even in this case the inverse pion 

u, 2 (Q 2 )) (where u a (Q 2 ) 



propagator, obtained from the l.h.s of eqs. (|4.5|) , is not equal to Z a (Q )(u) 



corresponds to the positive pion energies of isospin state a as given by eqs. (|4.5|) ) in an unambiguous way, 
since Z a has energy dependence determined only on-shell. Nevertheless, this last effect is suppressed in 
the non-relativistic limit. Indeed, C2 + c 3 is about one order of magnitude larger than g^/8mj^. 

Let us now introduce the couplings / t ° (temporal) and f® (spatial) as the value of the matrix elements 
(fi| Aq\-k°) /iQo and (£l\A"-\ir°) /iQi at threshold, i.e. for vanishing three-momentum. In the same way we 
also define the couplings / t + and / s + as the matrix elements (0| J^ \tt + ) /i\^2Qo and (fi| J^ 2 \ ir + ) /i\/2Qi 
at threshold and analogously for the tt~ with the couplings ff and /~, that can be obtained from the 
former just by exchanging p — ► —p. 

In the isospin limit, m = p = 0, the weak couplings for all the pions are equal and given by: 



ft 

fs 



2p 



P { 1 + 72 ( ° 2 + C3 



8m N J J 

)}• 



(4.25) 



see also ref. |T^]. From these equations we see that the spatial coupling, f s , is suppressed with respect to 
the temporal one by a factor: 

~ 2 _ n 4/5c 2 



1 



P 



(4.26) 



since C2 is positive as already discussed. Remember that v, given in eq. (|4~7|), corresponds to in-medium 
pion velocity in symmetric nuclear matter in the chiral limit. Taking into account the values for the 



low-energy constants given in sec. 4.4, we can write for ft and f s : 



U { 1 



Po 



-(0.26 ±0.04) 



}■ 
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Po 



-(1.23 ±0.07) 



}■ 



(4.27) 



with po the nuclear saturation density. We see that the reduction with density is sizeable for ft, being 
1/4 f n at po, and dramatic for f s , which vanishes already for p ~ 0.8po- We remark that in the vacuum, 
the study of the left-right current correlator leads to the conclusion that spontaneous chiral symmetry 
breaking is unambiguously signaled by the non-vanishing of the pion decay constant f n . By considering 
the coupling of the pion to the axial- vector charges Q A = J (fx Aq (x) , the generators of the axial chiral 
transformations at the quantum level, one can see that it is the temporal component ft that plays the 



role of f n in the medium. Indeed, a straightforward calculation taken into account eqs. (4.24) shows: 



(7r a (p)\Q k A \n) = / dxe- l » x {7r a \A k (0)\Q) = (27r) 3 5(p)iQoft 



(4.28) 



Mathematically, one can obtain a more meaningful result than eq. ( [4.28 ), that involves the product 
8(jp)iQo which in the chiral limit needs some care, by considering the wave packet: 



dp 



(2vr) 3 2|p 



■/(p 2 )k a ( P )) 



(4.29) 



with /(0) = constant as discussed in ref. |4q| . Thus, as long as ft ^ the vacuum is not invariant 
under axial chiral transformations and one has for symmetric nuclear matter the same chiral symmetry 
breaking pattern as in vacuum, namely, SU(2)l x SU(2)r — > SU(2)y. Of course, the fact that f s / f t in 



matter (and some consequences thereof) was already discussed in the literature, see e.g. [12,13,17 
In addition, it should also be interesting to perform the in-medium analysis of the left-right current 
correlator as in the vacuum. 

One can give a more microscopic picture for the large quenching of the spatial components of the 
axial-vector current by considering the resonance saturation of the 0(p 2 ) CHPT counterterms, the c's, 
making use of ref. [^J. In this way we also connect with earlier many-body approaches [48, 49), which 
stress the importance of the isobar A (1232) in the related quenching of Gamow- Teller transitions in 
finite nuclei. In fact we agree with this observation since the A contributes about 80% to the very large 
number C2 — C3 = 7.9 GeV -1 which controls the in-medium quenching of f s , eq. ( |4.27| ), while the other 

mm* 



main contribution, around a 15%, is due to the crossed exchange of a heavy scalar resonance 
On the contrary, the contribution of the A to ft vanishes because then the sum C2 + C3 appears. In this 
case the main contributions arise equally from the TV* (1440) and the crossed scalar exchange and amount 
to a much smaller quantity than the difference C2 



C3 as shown in eq. ( 4.27 ). 

Now, let us take into account eqs. ( [4.25 ) and the resulting expression, from eqs. ( 4.22j ), of the in-medium 
pion mass in the p— isospin limit, just defined above: 



Ml = Mi 



1 + ^ 
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2ci -C2-C3 + 



8m n 



)}■ 



(4.30) 



which, together with eq. ( |3.4[ ) for the p— quark condensate (£l\uu + dd\Q), establishes that in-medium 
corrections up to the next-to-leading order for the case of symmetric nuclear matter do not destroy the 
validity of the Gell-Mann-Oakes-Renner relation (GMOR), which is only exact at lowest order: 



ft Ml 



-m(ft\uu + dd\Q) + 5 V 



(4.31) 



#9 Although in ref. [g3| this contribution was mimicked by the exchange of a light a, in ref. [^0| it was pointed out that 
when chiral symmetric Lagrangians are used and unitarized in harmony with CHPT to describe nN scattering, then the 
resulting crossed scalar contribution comes out to be heavy. 



19 



where here 5 vac is a vacuum CHPT correction at 0(p 4 ). The stability of the GMOR relation under 
in- medium corrections as well as the fact that it is the temporal coupling ft, and not f s , that is the 
one involved in the GMOR relation, has been previously reported in refs. 12] , 13| , within the mean field 
approximation, and in ref. |HJ ], in the framework of QCD sum rules. 

We now consider the general case with p ^ and turn our attention to the pseudoscalar isovector 
quark current 

P i = qhs^q , (4.32) 

which can be obtained following similar steps as in the case of the axial-vector currents by differentiating 
the generating functional with respect to p l (x) with p(x) = X]iP* r *> taking finally the limit of vanishing 
external sources. We are interested in evaluating (Q|P*|7r fc ). Since the pseudoscalar source counts as 
0(p 2 ) it first appears in A^> and one has to evaluate an analogous diagram to that of fig.^ with the 
wavy line indicating now a pseudoscalar source. It is then straightforward to obtain: 

I r V 8m N Qq 

2 \<r2 



with the vacuum pseudoscalar pion coupling calculated at 0(p 4 ) Q. 

This leads naturally to the study of the Ward identity relating the divergence of the axial- vector quark 
current with the pseudoscalar ones, 

<9 M (qjvl^q) = mqij 5 T l q+ 5 l3 mqi-y 5 q . (4.34) 

Setting fh = and evaluating this expression at threshold one obtains a relation between f® ,+ '~ , M 2 0i+t _ 
and G%, and G~ . The latter are defined in terms of the matrix elements eq. ( |4.33 ) at threshold, 



when the factor \pl is removed in the charged matrix elements. Considering first the 7r° case and from 
eqs. (p>4D and (|4T33|) , it follows: 



= UMl 1 1 + ^ U C1 -c 2 -c 3 + \=mG n , (4.35) 

since of course in the vacuum f v M 2 = mG^ holds ||. Then 

f t Ml« = mG n , (4.36) 



see also ref. [12]. An analogous relation for the charged pions can be also obtained by considering the 
divergence of the axial current J^ 2 together with the pseudoscalar current P 12 . One can readily check 
that 

f t + M 2 + = mG n , 

ffM*. = mG^ . (4.37) 
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It is worthwhile to stress that in eqs. ( 4.36| ) and (4.37) we have not imposed that p = but just fh = as 
is necessary in order that the Ward identity of eq. ( 4.34 ) holds in terms of only the isovector pseudoscalar 
currents. 

We now consider the general case with To study eq. (4.34), it is necessary to consider as well 

the isoscalar pseudoscalar current: 



P = qij^q = UIJ5U + dij^d 



(4.38) 



that can be obtained from the action J dyC^ in the same way as before by differentiating with respect 
to p°(x) such that p(x) = ^2iP 1 {x)T' 1 -\-p°I with i" the 2x2 identity matrix. It is straightforward to check 
that 



^-3 



8Bc 5 



p-S 



8B 2 



t 7 fh = 8 h3 G 



(4.39) 



/ r f 

where we have explicitly shown the vacuum contribution proportional to i 7 from ref. 0]. Proceeding 
along analogous lines as those to derive eq. ( [4.36 ), when considering the product / ( M^ one has to add 
to mGj,- the term: 

raG:+myC 5 , (4.40) 

where the first term comes from the vacuum equality f n M 2 = fhG^ + fhG* and the second from the 
perturbative expression of the ir° mass, eq. ( 4.22j ). Collecting both terms and comparing with eq. flOp 
one finally has: 

/ t °M^o = mG^ + fhG* . (4.41) 

Note that the relations for the charged pion quantities in eqs. ( [4.37 ) do not get any contribution from 
fh 7^ at the next-to-leading order considered here. 

We have also checked that the Ward identity eq. (4.34) is satisfied as well for non-zero three-momenta, 
that is, above threshold. This constitutes a good check for all our former expressions, namely: the 
spectral relation between energy and three-momentum, eqs. ([4.5|), the matrix elements of the temporal 
and spatial components of the axial- vector currents with one pion field, eqs. ( 4.24j ), and the pseudoscalar 
ones eqs. ( gjg ) and flpgp . 



5 Three-point Green functions 

In this section we will discuss the in-medium contributions, up to 0(p 5 ), to the coupling of pions to a 
vector source. As will be shown below, the conservation of the electromagnetic current is nontrivial in this 
context. Moreover, the p— corrections at threshold are of the order of 100% of the vacuum ones already 
at p = po due to a large counterterm contribution proportional to C2. Nevertheless, for higher energies, 
the p— corrections become considerably smaller compared to the vacuum ones due to the dominant role 
of the counterterm £ G of SU(2) L xSU(2) fi CHPT g. While c 2 is saturated by the A(1232) @ resonance, 
£q is saturated by the p(770) ||. Furthermore, we will show that the surrounding nuclear medium does 
not alter the anomalous 7r° — > 77 decay amplitude up to 0(p 5 ). We will also briefly discuss the vanishing 
of three pion scattering in the medium. 



5.1 Pion coupling to a photon source 

Let us study the process 7* — > 7r l (q)7r J 1 (— q). In the effective field theory the photon field A^(x) is included 
via the external vector source v fl (x) = eA fJi (x)Q, with eQ = ediag(l,0) the nucleon charge matrix. The 
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Figure 6: In- medium diagrams with two pions and one photon up to 0(p 5 ). The analogous diagrams with more 
than one Fermi-sea insertion are not shown, as they do not contribute here. 

in-medium contributions are depicted in fig. ||. In those diagrams where two and more insertions 
are indicated, only one of them has to be considered in the calculation up to 0(p 5 ). Indeed the only 
non- vanishing contribution with an vertex is that of fig.^a, the other ones are at least 0(p e ). The 
calculation is straightforward and for on-shell pions gives: 

(n+n-\V \n) = e j--^ + ^p + W+ (q)-u,_(q)} , (5.1) 
(vr vr \V |fi> = -2eq F y (q ) jl - - _ j2 + _ j2 + -^jr^j , (5-2) 

where F\/(q 2 ) is the vacuum pion vector form factor up to 0(p 4 ) given in ref. ||. Note that there is no 
coupling to a 7r°7r° state as expected from Bose-Einstein statistics. The term w+(q) — w_(q) in eq. ( ^T| ) 
arises because of the difference in the energy of the ir + and n~ in asymmetric nuclear matter as given 
by eq. (4.5). Let us stress that the conservation of the electric current, c^V^ = 0, is nontrivial for the 



in- medium case, as one can find out by studying the origin of the different terms in eq. ( |5.l|) . The first 
one comes from fig. ^a, the second from figs.|6|f, g and h while the last two terms originate from from 
graphs with only two pion lines, without photons, see eqs. (|4.5|). Despite of that, all of them cancel each 
other and the temporal derivative is zero as required to guarantee the conservation of the electric current. 
Indeed if one calculates the matrix element (7T + 7r — IcLV^lO) from eq. (5.1) and (|5,2j ), one obtains: 



(n + TT-\d^\n) oc 2^q)(--^ + ^p) +u; + (q) 2 -^(q) 2 



P ) \ P ^ /Mq), 

= . (5.3) 

The in-medium contributions to the spatial components of the vector current, eq. (|5.2|), come from figs.^a 
for the first density-dependent term inside the curly brackets, from |||b, c for the second one, and from 
fig.||f, g, h for the last two terms. Since it does not involve the non-relativistic suppression factor l/mjy, 
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the dominant term is the one with the rather large C2 counterterm which is also enhanced by a factor 



of 2 as compared, e.g., to eq. (4/7). The third term in eq. ( |5.2D , figs.yb and c, involves only the proton 
density since only protons couple to the photon at this order. The contribution from figs.|6]d and e are 
canceled by those of the wave function renormalization. The last term comes from fig. ||g and this is why 
it only involves the neutron density. The term before the last one is a combination of figs.^f and h, with 
only proton density dependence, and that of fig. pig. 
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Figure 7: In-medium pion vector form factor in symmetric nuclear matter as function of the modulus of the 
pion 3- momentum (= |q|). The continuum line is the vacuum 0(p 4 ) result from ref. ||. The dashed, dotted 
and dashed-dotted lines are the in-medium results up to 0(p 5 ) and correspond to kp — k F /2, 3fc^/4 and k F , 
respectively. Note that at q V7r — the photon virtuality is four times the pion mass squared and that the vacuum 
form factor is one at virtuality zero. 

In fig.0 we show the in- medium results to the spatial vector coupling in symmetric nuclear matter for 
three values of kp: k F /2, 3k F /A and k F . As stated above, the modifications with density increase 
rather fast due to the large C2 counterterm. Note, however, that the relativistic corrections become more 
important for higher three-momenta. In fact, the difference between the vacuum 0(p 4: ) result, solid line, 
with respect to any of the other lines, gives the in-medium contributions. They rapidly increase as soon 
as we approach k F because of the cubic dependence on kF over a quantity of order k F . The lowest 



order 0(p 2 ) result is just 1. Taking the resonance saturation of the C2 into account, see ref. [23|, one 
can see that this counterterm just reflects the impact of the A(1232) in low energy ttN scattering. Thus 
the term proportional to C2 in eq. (^^), stemming from fig.^i, establishes a large A-hole contribution. 
Nevertheless, for larger values of the three- momentum |q|, the in- medium corrections are relatively less 
important than the vacuum ones, due to the large £q CHPT counterterm that signals the well-known 
and very important p(770) contribution to the pion vector form factor in vacuum 0. Finally note that 
the form factor in nuclear matter does not go to one at zero virtuality, because this background carries 
charge. 
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5.2 Anomalous sector 



Since the effective Lagrangian C nn contains terms with odd number of pions including a 75 matrix, one 
might expect the in- medium appearance of processes violating intrinsic parity, as, e.g., the scattering of 
three pions or contributions to the 7r° — ► 77 decay amplitude. We will show that, up to 0(p 5 ), there 
are no in-medium corrections to the ir° — ► 77 amplitude and that the three-pion scattering vanishes in 
unpolarized nuclear matter at the same order as well. 

The set of diagrams to be considered are shown in fig. [1. The figure does not contain any diagram with 
only one A vertex, as such a contribution is trivially zero. Namely, when A = A^\ it vanishes because 
then we have one 75 matrix together with, at most, two other 7 matrices, such that the Dirac trace is 
zero. When A = A^ 2 \ the result is zero, simply for the reason that there are no terms in A^ with two 
photons and one pion. As indicated in the diagrams of fig.|§], each diagram can contain an operator 
attached to a photon line. In fact, one could also have diagrams with just two photon lines stemming 
from the same A^ vertex; but they vanish since [v fl (x),v l/ (x)\ = 0. 




a) b) c) d) e) 



Figure 8: In-medium diagrams with two photon and one pion lines up to C(p 5 ). The analogous diagrams with 
more than one Fermi-sea insertion or with only one A vertex are not shown. Every diagram vanishes independently 
up to 0(p 5 ). 

Let us show that all the diagrams of fig.|8| vanish. Consider first the diagrams fig.||a and b with only 
two leading A^ operators. It is easy to prove that the flavor trace vanishes for each of them. For this 
case T^(x) = — iv^x) and A At (x) = — [(p(x) , v^(x)]/2. As a result the flavor traces in figs. |8|a,b vanish 
since: 

n(p)(r M (x)A„(i/)) -> (Q[0,Q])=O, 
n(p)(r 3 T^x)A u (y)) -> (r 3 Q[0, Q]) = , (5.4) 

where Q is the charge matrix. A change of the order of T^(x) and A u (y) in the equations above does not 
modify the results since [73, Q] = 0. For the same diagrams, but now with one A^ operator replacing 
one A^ operator, there is only a contribution when the other A^ = —igAl^ls^-n-i with A^ as given 
above. p 1 ^! For the terms in A^ without gamma matrices the result vanishes, after taking the Dirac 
trace, because of the 75 accompanying A^ in A^ and the fact that only three extra gamma matrices are 
available. The contribution proportional to C4 in A^ vanishes because there are no terms in involving 
vector sources without pion fields. Hence two or more pions would arise and this is not allowed. The 
terms with cq and C7 involve the which for our present purposes simplifies to 2d li v v — 2d u v fl . Thus 
they only contain the charge matrix Q in flavor space. In summary, when is replaced by F^ u , the 
same reasoning as previously given for the corresponding diagrams with only A^ operators holds also 
for this case. The diagrams fig. ||c, d and e deserve special attention, although in the end each of them 

* 10 The simplest way to realize this is by considering the ip— counting introduced in ref. Q where <j> ~ ip, a M ~ tp, p ~ ip, 
Vfj, ~ ip 2 and s ~ ip 2 . In this counting, the operator always gives rise to terms with odd powers ol ip and the F M operator 
as well as to even power terms. 
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also vanishes independently at 0(p 5 ) because of the Dirac trace after symmetrization with respect to 
the photon lines. Furthermore it is easy to see that the diagrams fig.||c, d, and e with one vertex, 
attached to a v^{x) external source, vanish as well. Since these diagrams are at least 0(p 5 ), one can 
replace the momentum pj by p inside the expression j/j +m,N, which enters in the numerator of the baryon 
propagator, because pj = p + Qj and the difference will be 0(p e ). Furthermore, for the same reason, we 
can set j/ equal to jom^- Taking also into account that the four-momentum Q of the ir° is just (M°, 0) it 
is straightforward to see that in all the Dirac traces one has the product 7o(I + 7o)7s(-f + 7o) = which 
implies the vanishing of the 0(p 5 ) contributions from such diagrams. 

As a result of this discussion, taking also into account that there are no contributions from diagrams 
with two or three Fermi-sea insertions because of the reasons discussed in sec. 4.1, the final in-medium 
71-0 — * 77 decay amplitude T(ir° — ► 77), up to 0(p 5 ), just corresponds to the pure vacuum one given by 
the Wess-Zumino-Witten term in the effective chiral action [j52]| , counted as 0(p A ) in CHPT. Note also 
that the corrections to the vacuum phase space induced by the change of the 7r° mass in the medium 
are at least 0(p 9 ), of the same order as the neglected contributions to the decay resulting from 0(p e ) 
contributions to T(-7r° — ► 77). Thus up the accuracy we are considering in this paper, there is no 
modification of the tt° —* 77 width due to the nuclear medium. 

Finally, let us consider the 3tt scattering, which is not allowed in vacuum, since the anomalous WZW 
terms only give rise to processes with five or more Goldstone bosons. It turns out that the 3tt scattering 
does not occur in nuclear matter, at 0(p 5 ), either ^ TT . The set of diagrams to be evaluated is depicted 
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Figure 9: In-medium diagrams with 3ir lines up to (D(p 5 ). Each of them vanishes independently in asymmetric 
nuclear matter. The analogous diagrams with more than one Fermi-sea insertion are not shown. 

in fig. [| The diagram of fig. ||a vanishes trivially because of the Dirac trace since we have one 75 and 
at most two gamma matrices while at least four gamma matrices are necessary. The latter structure is 
possible for diagrams fig.|9]b,c,d; but in this case one is left with the totally antisymmetric tensor in four 
dimensions, e^ vpa , which requires four independent momenta. However, we have only three of them in 
any diagram: one is the in-medium running on-shell four-momentum p together with two four-momenta 
from the pions (the third one is given by energy-momentum conservation). Thus these contributions 
vanish. Note that only diagrams with just one Fermi-sea insertion are shown in fig.||, since the proof 
that each diagram vanishes independently is only based on the Dirac structure of the diagram which is 
the same independently of whether a free baryon propagator is replaced by a Fermi-sea insertion or not. 



6 Four-point Green functions: Pion scattering 

The in- medium An scattering contributions begin to appear at 0(p 3 ) because an in-going pion four- 
momenta can be combined with an out-going one, such that Q® ~ for one of the propagators. The 

* n These will be subprocesses in the discussion of the more realistic 4tv scattering in the next section. 
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corresponding diagrams are shown in fig. [T(| One can see the 0(p 3 ) behavior just by applying eq. $L%, 
case b). If all the four external pion legs had Q° ~ k F /2rriN, then the chiral power counting would start 
at v = 1 with V p = 1. However, since the pion legs are of course on-shell with Q° > M n , fig.|lCi|a starts at 
0(p 3 ). Indeed, the analysis of all the possible diagrams according to eq. (|2.9[)-b) up to 0(p 3 ) shows that 
the leading contribution corresponding to the diagrams shown in fig.[l(] appears at v = 3. 




a) b) c) d) 



Figure 10: In-medium diagrams contributing to An scattering at 0(p 3 ). All the vertices are operators. 



Only one Fermi-sea insertion has to be considered in fig. 10 
four of them, because there is always one four-momentum Q° 



makes it impossible to satisfy pj 



in 



N 



There is no way to include three or 
M v running along one insertion that 
for this thick line, as discussed in sec. |4.1|. However, it is in 



principle possible to have two Fermi-sea insertions. Nevertheless, as shown in sec. |4.1| , in this case there 
is a cancellation, below the Fermi-momentum, between the diagrams with two Fermi-sea insertions and 
the imaginary parts of those with only one Fermi-sea insertion. It is easy to see that for Q° = both 



contributions cancel each other: setting pj = p + Q and then imposing p - 

Q 2 + 2pQ = , 



m 



N 



, one can easily see that: 
(6.1) 



■m 



AT" 



On the other hand, because of the previous cancellations, |p + Q| > kp. By 

> kp that is not allowed. 



since Q° = and p 2 

squaring this condition and imposing eq. (|6.1| ) , one concludes p 2 

In fig. [ll| we show some 0(p i ) diagrams that are two connected 3ir scattering processes, which vanish, 
as already discussed in the previous section, because of the appearance of the totally antisymmetric 
tensor e^ uap . Nevertheless the vanishing of the diagrams in fig. 11 is due to the pseudoscalar character 
of the exchanged pion. If instead there were short range NN forces, these diagrams would not only 
have survived, but would be enhanced by the appearance of nucleon mass factors from each bubble, as 



discussed in sec. 4.2. 






a) 



b) 



c) 



Figure 11: In medium-pion diagrams contributing at 0(p 4 ). Each of them vanishes independently. More details 
are given in the text. 

In the following we will restrict ourselves to the case of symmetric nuclear matter p = 0. The results 
presented here have been calculated in two ways: using the relativistic formalism from the beginning 
and performing the chiral expansion up to C(p 4 ) afterwards or, from the beginning, making use of those 
simplifications that are allowed order by order. The results are the same under both ways. Because of the 
loss of crossing symmetry due to the absence of Lorentz covariance, we consider separately the scattering 
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amplitudes: 



T 2 = TT + (q 1 )TT + (q 2 ) -> TT + (q 3 )TT + (q 4 ) , 

T 3 = TT + ( qi )<ir-(q 2 ) - 7r + ( g3 )7r-(9 4 ) , (6-2) 

where % is the four- momentum of each pion with qi + q2 = q3 + q4 = 0. We further define the three- 
momenta Q = qi — q 3 and R = qi + q 3 together with the usual Mandelstam variables s = (qi + q 2 ) 2 = 
(Q3 + <?4) 2 , t = (qi - <? 3 ) 2 = (q 2 - qi) 2 = -Q 2 and u = (qi - g 4 ) 2 = (q 2 - q 3 ) 2 = -R 2 . It is worthwhile to 
mention that the 0(p 4 ) corrections to the leading 0(p 3 ) diagrams in fig.[l(] vanish. Note as well that the 
in- medium corrections to the pion mass in the symmetric nuclear matter, eq. ( |4.22 ), as well as the wave 



function renormalization, eq. ( D.14j ), appear first at 0(p 5 ). Then up to 0(p 5 ) and in the p— isospin limit 



m = p = 0, one has the final expressions: 

rV^-^O = ^ + ^^[J(0,Q)+T(0,R)] , 

U.V^nV) = ^^-^[X(0,Q) + X(0,R)] + ^{(t + 2. 2 -2M 2 )T(0,Q) 

+ ( U + 2u; 2 -2M 2 )T(0,R)} + ^^{^-2( W 2 -M 2 ) 2 }[J(0,Q)+J(0,R)] , 

liCir+^-^O = ^^ + ^X(0,Q)-^{t + 2. 2 -2M 2 }X(0,Q) 

+ ^^(- 2 -M 2 ) 2 J(0,Q), (6.3) 



with lo the energy of any of the pions and the function 2(Q , Q) as defined in eq. ( |C.2|) . Notice the presence 
of the 2itin factor accompanying the 1(0, Q) and 1(0, R) functions that, as discussed in sec. 4.3, implies 
a dramatic reduction of the chiral scale A to about 170 MeV. Nevertheless, one should keep in mind that 
this scale only refers to the Fermi-momentum kp and the three- momenta Q and R, but not to the energy 
oj of the pions which is bounded from below by twice the pion mass. The contributions in eqs. ( |6.3; ) 
proportional to g\ and g\ come from fig. ICa and fig.[lO|b,c, respectively. Those stemming from fig. |l0|d 



originate from the Weinberg- Tomozawa term and do not contain any g& factor. 

From the previous expressions one can easily work out the in-medium amplitudes with well defined 
isospin (I) and angular momentum. These are shown in figs. [l^ where, from top to bottom, the first 
panel corresponds to the S-wave 1=0 amplitude, the second one to the P-wave 1=1 partial wave and the 
third panel to the S-wave 1=2 amplitude, for three values of the Fermi-momentum fej? as indicated in the 
figure. The thick solid curves correspond to the lowest order tttt CHPT partial waves while the thin solid 
ones to the pure vacuum 0(p 4 ) tttt amplitudes, ref. 0. One can see strong in-medium corrections close to 
threshold for the 1=0 and 2 S-waves induced by the diagram of fig. [Hjd. This dia gram only contains the 
Weinberg- Tomozawa term proportional to T^, eq. dB.3[), which does not involve a pion three-momenta 
as is the case for the term proportional to in eq. ( [B.3D . This diagram has so far not been considered 
in the works (see [p3-55|) that are involved in the physics of in- medium pion scattering, triggered by 



the experimental data of the CHAOS Collaboration ||56j. The latter reported originally a dramatic 
enhancement of the (7r + , 7r + 7r~) reaction close to the 2ir threshold, although opposite conclusions have 



recently been obtained by the Crystal Ball Collaboration [57] from the study of the reaction (n , n ir ). It 



is nowadays believed that when one properly accounts for the limited acceptance of the CHAOS detector, 
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Figure 12: From top to bottom: S-wave 1=0, P-wave 1=1 and S-wave 1=2 7T7t partial wave amplitudes. The thick 
full lines are the leading 0(p 2 ) tttt scattering amplitudes while the thin solid ones are the pure vacuum 0(p 4 ) 
results from ref. ||. The dashed, dotted and dashed-dotted lines are the 0(p 3 ) in-medium partial wave amplitudes 
for kp — 100, 200 and 300 MeV, respectively. 



these results are mutually consistent As one can see, the diagram fig.[To|d is indeed dominant for small 
pion three-momenta, as far as this perturbative analysis holds, since it is not kinematically suppressed as 
the other contributions with qa dependence because of the P-wave coupling of the pions to the nucleons. 
At this point it is also worthwhile to note that we have shown explicitly, eq. (|6.3D , that the in-medium 
contribution in symmetric nuclear matter to tttt scattering from the leading pion-nucleon Lagrangian, 
C^lf, does not vanish as stated in ref. |5f|. On the other hand, our results for the S-waves, valid for 
nuclear matter, cannot be taken literally for finite nuclei when Q° = 0, Q = since theI(Q ,Q) function 
should then vanish due to the finite excitation energy of a few MeV of the nuclear energy levels. This 
effect is pointed out in ref. |59[. This would imply that in finite nuclei the S-wave 1=0 and 1=2 amplitudes 



at threshold do not obtain any contribution from fig. lCd, although, as soon as the pion three-momentum 
is different from zero, these large corrections rapidly build up. However, as energy increases, the diagrams 
of fig. lCb and c, which were as well not considered in ref. [p3|-|55[, become more relevant, dominating 
the in- medium corrections to the P-wave 1=1 scattering amplitude. On the other hand, their influence 
is rather small for the S-waves. The term proportional to g\, fig. |lC|a, absent as well in refs. 53-55|, 

#12 See the last entry of ref. |59]. However, this issue is still controversial pa]. 
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is very relevant for the S-waves when the pion three-momentum is higher than 100 MeV, although its 
influence is much smaller for the P-wave. It is remarkable to point out that for the 1=0 and 1=1 cases 
the amplitudes change dramatically with density, while for the 1=2 case the S-wave scattering amplitude 
receives in-medium corrections smaller than those from vacuum 0(p 4 ) CHPT when |q| > 200 MeV. 
Moreover, the 1=0 partial waves changes even sign, becoming repulsive, while in the 1=1 case there is a 
strong enhancement as soon as one moves to higher three-momenta and Fermi-momentum. We remind 
the reader that in the present situation with Q® ~ the perturbative results are only valid for )zf and 
three- momenta up to around 170 MeV and that the dotted and dashed-dotted lines in fig. ^ are just 
shown to point out the huge in-medium corrections as soon as one goes beyond that limit, although 
in the end the actual size of the in-medium corrections depend on the specific reaction. Indeed, for the 
physics relevant to the CHAOS experiment the effective average density is around po/2 [^] corresponding 
to a kp ~ 210 MeV. Hence, one would expect large in-medium non-perturbative effects to take place. 
Nevertheless, the different models dealing with this problem should match at low densities and momenta 
with the 0(p 3 ) results presented in eq. ([T^) in order to take care of the requirements of chiral symmetry. 



7 Conclusions 



In this paper we have tackled the problem of establishing a chiral effective field theory in nuclear matter 
with explicit pion fields. We have made use of the results of ref. (l[ where, by integrating out the baryonic 
fields in the path integral representation of the generating functional, once the change of ground state from 
vacuum to nuclear matter is done, the in-medium chiral SU(2)xSU(2) Lagrangian is derived. We have 
briefly reviewed these results in appendix |A| as well as the perturbative case presenting the corresponding 
Feynman rules. After establishing the power counting rules, we have systematically studied several low- 
energy QCD Green functions up to next-to-leading order, (D(p 5 ), when the standard counting holds or by 
working out the leading in-medium contributions in the non-standard case. The novel results obtained 
here can be summarized as follows: 

(1) In contrast to previous works, which apply the mean-field approach or many-body calculations, 
the in-medium chiral counting is worked out in sec. 2.1, The counting scheme is dependent on 
the energy flowing into the nucleon lines. This leads one to consider the standard and the non- 
standard case, respectively, as summarized in table ^. In the former case, the chiral expansion of 
pion properties in the medium starts with terms at 0(p 4 ), and the next -to-leading order corrections 
appear at 0(p 5 ), quite different to the in-vacuum power counting. 





Standard counting Non-standard counting 


energy flow 
nucleon propagator 
counting index v 
breakdown scale A 


Q° ~ M w ~ 0(p) Q° ~ Q 2 /2m N ~ 0{p 2 ) 
D, 1 ~ 0(p~ l ) D l ~ 0(p- 2 ) 
eq. Q eq. Q 
VSirU ~ 0.7 GeV 6ir 2 f 2 /2m N ~ 0.27 GeV 



Table 2: Chiral counting for in-medium CHPT for the standard and 
the non-standard case. The pertinent breakdown scales are also given. 
Here, "energy flow" means the energy flowing into the nucleon lines 
mediated by pions or external sources (as defined in sec. [O]). 



(2) We have also established the relevant scales of the problem when restricting ourselves to the C wn and 
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C-xN Lagrangians, see sec. ^ or ref. [||. In the vacuum, the pertinent scale is A x ~ 1 GeV~ Airf^. In 
the medium, one has two new scales. These are: VW/tt — 0.7 GeV and Qtt 2 f% /2m^ ~ 0.27 GeV, 
in case that the standard or the non-standard counting rules apply, see also table |2|. We point out 
that in case of P-wave interactions, these scales are reduced by factors 1/gA and l/g\, respectively. 

(3) We have studied the quark condensates and re-derived, from the effective field theory point of view, 
known results in symmetric nuclear matter, and have further extended them to the non-symmetric 
case. 

(4) We have considered the propagation of pions in the medium obtaining the pion propagator up to 
0(p 5 ). In this way we have established that chiral symmetry can account for the observed shift 
of the mass of the negative pion in deeply bound pionic states in 207 Pb. Our numerical result 
AM n - = 18 ± 5MeV is compatible with the experimental number, AM W - = 23 - 27MeV §,0 
within errors. 

(5) The wave function renormalization of the pion fields corresponding to the calculated in-medium 
action J dxC has been established making use of first principles. 

(6) We have also studied the coupling of pions with axial- vector and pseudoscalar sources. In particular, 
it is shown that in- medium corrections up to 0(p 5 ) do not spoil the validity of the Gell-Mann- 
Oakes-Renner relation. We have also found a decrease with increasing density for both the quark 
condensates and the temporal component of the pion decay constant ft- Both effects seem to 
indicate a partial chiral symmetry restoration with increasing density. We remark again that a 
systematic study of the in-medium order parameters still has to be performed. A drastic quenching 
with density has also been obtained for the spatial component of the pion decay constant f s . To 
0(p 5 ) we have checked the QCD Ward identity relating both the temporal and spatial components 
of the axial- vector currents with the pseudoscalar ones and quark masses. 

(7) A rapid decrease with density of the coupling of a photon to two pions, particularly in the threshold 
region, has been found. The derived vector current amplitudes coupled to two pions fulfill the 
requirement of current conservation. Furthermore, we have established the absence of in-medium 
renormalization up to 0(p 5 ) of the anomalous 7r° — > 77 decay amplitude. 

(8) Finally, tttt scattering has been studied up to 0(p s ) since in this case the non-standard counting 



occurs. As explained in sees. 2J. and y, this implies that the in-medium corrections start at lower 
orders than in the standard case, here already at 0(p 3 ). In addition the scale, below which the 
perturbative expansion is applicable, decreases. As a result the in-medium corrections increase very 
rapidly with density and already at kp ~ 200 MeV, or at a density of just ~ 0.4poj they are 100% 



with respect to the lowest order CHPT results. The diagrams presented in figs. 10 are the leading 
ones in the chiral expansion and have not been considered so far in the literature. 



Future challenges are the inclusion of multi-nucleon contact interactions, which are enhanced because 
of the largeness of the S-wave scattering lengths related to the presence of shallow NN bound states, as well 
as the simultaneous calculation of pion loops necessary for determination of the full C(p 6 ) contributions. 
Furthermore, the possibility of some non-perturbative scheme that allows for the recovery of the scale 
V^vr/^, even in the case of the non-standard counting or when the multi-nucleon local interactions are 
included, should be pursued. 
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A Generating functional 

Our starting point is the generating functional Z(v, a, s,p) in the presence of external vector v, axial- 
vector a, scalar s and pseudoscalar p sources (included in C nn and C^n) which is obtained after integrating 
out the baryonic fields, as derived in ref. Q (for all details on the derivation, we refer to that paper). It 
is given by 

e iZ[v,a,s,p] = J ^ exp | -y dx c ^ 

1 f dp f dq f dxdx > dyd y> e ip(x-y) e -i<i(x'-y') (A.l) 



2 J (2ir) 3 2E(p) J (2Tr) 3 2E(q) , 
x Tr\A[l4 - D^A]- 1 ]^ {({ +m N )n{q)A[I A - D^A}' 1 ] 

(y',y)(P ' + m N)n(p) S j + ...| , 

where the trace, indicated by Tr, is over the spinor and flavor indices. From this result we can readily 
read off the new effective chiral Lagrangian density C n7T in the presence of nucleonic densities just by 
equating the expression between curly brackets to i J dxC nn . Furthermore, the diagonal flavor matrix 



n(p) in eq. ( A.l ) is: 

' 0(k ( P " IP 

0(k™-\p\) J \ n{p) 




with I2 the 2x2 unity matrix, 0(x) the Heaviside step function, T3 the usual Pauli matrix, T3 = diag(l, — 1), 
h(p) = (rt(p)i + 71(^)2) /2 and n(p) = (n(p)i— n{p)2) /2. We are considering the general case of asymmetric 
nuclear matter with two Fermi-seas of protons and neutrons with densities fP = (kp^) 3 /3ir 2 and p n = 
(k^) 3 /3tt 2 , respectively, with kp and k^ the corresponding Fermi momenta. In eq. ( |A.1| ) the nucleon 



energy is given by E(p) = yp 2 + m? N , where tun = {n^n + m p )/2 — 939 MeV, since the difference 
between the proton and neutron masses is 0(p 2 ) giving rise to in-medium contributions of higher order 
than the next-to-leading order ones, see section |2.1| . For the same reason, we will also use the average 
nucleon mass, mjv, in the baryon propagators. 

The ellipses in eq. (A.l) indicate terms with a higher number of three- momentum integrals over the 



Fermi-seas coming from a logarithmic expansion of the vacuum interaction operator T = — iA[I± — D^" 1 ] -1 
[Q]. Notice that the Fermi-sea states are on-shell, p° = E(p). Indeed, for n > 1 Fermi-sea insertions one 
picks up a factor (— l) n+1 /n, while every operator T = —iA[l4 — Dq 1 ] -1 , containing both pion legs and 
external sources, implies a phase of — i . The operator A is defined as follows: 

/ dx Cjtn = / dxip(x)D(x)ijj(x) = I dx ^(x)Dq(x)^(x) — / dxdyip(x)A(x,y)ip(y) , (A. 3) 
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where Dq{x) = i^f^d^ — is the free Dirac operator and the symbol Mjy refers to the diagonal matrix 
Mn = diag(m p ,m n ) of physical proton and neutron masses, respectively. Hence, standard in-medium 
generalized vertices from eq. ( A.l) can be represented as in fig. p] where each connected diagram represents 



one of such vertices and the vacuum non-local vertices T are connected through the exchange of on-shell 
Fermi-sea states. The latter are represented by the thick solid lines in that and all following diagrammatic 
figures. Furthermore, V results from the iteration of the — iA operator with intermediate free baryon 
propagators iD^ 1 as shown in fig.|2|. Consequently, any diagram with medium contributions will be a 
set of V p > 1 generalized in-medium vertices. For each of them one has rij > 1 Fermi-sea insertions, 
1 < J ' < V p , rrij > free baryon-propagators and of m,- + nj vertices — iA. The resulting Feynman rules 
(as derived in Q) read: 

• First include the global sign (—1)^ because of the fermionic closed loop attached to each generalized 
vertex, and the combinatoric sign factor 

from the logarithmic expansion giving rise to the different in-medium vertices. 

• Then, following the vertex in the opposite sense to that of the fermionic arrows, write for each 
Fermi-sea insertion an integral 



(2vr) 3 2£(p) 

with p° = E(p) and for every vacuum baryon propagator with free momentum p write 

. f dp y+rwjy 

V (27r)V-™^ + *e ' 1 j 

For a vertex in momentum space write a term 

-M(2vr) 4 , (A.7) 
keeping in mind the energy-momentum conservation at each vertex. 
Then for each pion internal line write the pion vacuum propagator 

dQ 1 (A.8) 



(2vr) 4 q 2 -M%+ ie 

• Finally, sum over the Dirac and spin indices of the fermions. 

This defines explicitly the Feynman rules in momentum space in order to obtain i(27r) 4 times the desired 
connected graph accompanied by the global delta function of energy-momentum conservation. It is 
important to keep in mind that a generalized vertex behaves as a standard local one of quantum field 
theory with respect to the determination of the numerical factors accompanying a given diagram under 
the exchange of pions. This can be seen by simply applying standard perturbative techniques in path 
integrals to the action given between brackets in eq. (A.l). It is also worthwhile to note that one can 



take eq. ( A.l ) directly in configuration space, see sections 



and O and the appendix 
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B Chiral Lagrangians and the interaction operators and 

To evaluate the next-to-leading order medium corrections to the Green functions from the irN Lagrangian 
£-wN = ^D(x)tp, as discussed in detail in sec. one has to consider the 0(p) C^L and 0(p 2 ) C^ N 
Lagrangians. For relativistic spin-1/2 fields chirally coupled to pions and external sources, the lowest 
order effective Lagrangian reads p2| 



4* = ^ m N I 2 + i-fT* + i 9 A 7^75 A J V , (B.l) 



with A M = \v)V v) in terms of the covariant derivative V ^U{x) = d^U{x) — i{v^{x) + a /1 (x))U(x) + 
iU{x){v jl {x) — a^{x)). is the chiral connection, = ~ [u^c^-u] — -1^(1^ + a^)u — 4ii(i^ — a^v) . The 
Goldstone bosons (the pions) are collected in the 2x2 unitary matrix u = exp(ii^/2/), U = u 2 and <p is 
given by: 



V2tt~ -tt° 



(B.2) 



On the other hand, the constants mj^, 5a and / refer to the mass, axial coupling of the nucleon and 
weak pion decay constant in the SU(2) chiral limit. Taking into account the definition of the A operator, 
A(x) = Dq(x) — D(x),^ T3 with Dq(x) = i^d^ — M^r, we obtain from eq. ( |B.1| ) for the first term in the 



chiral expansion of the operator A, 

A (l) = _^ F ^ _ i °g A . (B.3) 

The difference m N I 2 — Mjy is 0(p 2 ) j2^| and will be included in A^ 2 \ Although strictly speaking 
there is no explicit symmetry breaking at leading order in the effective pion-nucleon Lagrangian, this 
identification will turn out to be convenient at later stages. 

To obtain the complete expression of the (second order) A^ 2 ' operator, we have to take into account 
the Lagrangian C^ N (we use here the notation of ref. Q]). The operator A^ can be obtained by just 
removing the ijj and ip fields from the Lagrangian and by changing the overall sign, 

A< 2 > = m N I 2 -M N - Cl ( X+ ) + -^(u^D^D" - ^ + ^-f-fKM - c 5 ±+ 

2m N 1 4 

- -^tY^-^tY^), (B.4) 

8 m^v 8 m 

where the 0(p 2 ) low-energy constants Cj are finite and where we have also added the term m^l 2 — Mjv 
as discussed above. Furthermore, we have = + T^, = 2iA^, \+ = u*xu + ux'u with X ( x ) = 
2Bq(s(x) +ip(x)). Here, BqS^ = — (0|(f q J |0)/,/' 2 measures the strength of the symmetry breaking in the 

chiral limit. The quark niciss matrix J\A — d.iag(?n w , nid) is included in the scalar source s(x) = M H . 

The trace in flavor space of an arbitrary matrix O is denoted by (0). In this way, the traceless matrix 
X+ is defined as x + = X+ ~ Finally, F+ = u^F^u + uF^u f where F^ v = d^ v - dj^ - i[tp,l v ], 

Fjw = d^r u — d u r^ — i[r^, r v \ with = — and = v fJL + a^. For more details see e.g. ref. p4) . We 
briefly collect the chiral power (or chiral dimension) of the basic ingredients that appear repeatedly in 

#13 Although the operator A is in general non-local, A(x,y), we omit here the second argument since it just involves the 
delta function, 8(x — y), or derivatives thereof. 



33 



the chiral expansion of the operator A. Namely, d^u, and are ~ 0(p) while s, p, x and F^ R are 
~ 0(p 2 ). On the other hand, u and d^ij) ~ 0(1) (although diip ~ 0(p) for the spatial components i). 



We point out that the first two terms in eq. ( |B.4| ) appear because the free Dirac operator Dq is defined in 
terms of the physical nucleon masses (and not the mass in the chiral limit). It is straightforward to show 
that up to 0(p 2 ) in A@\ the difference I2 — Mn is canceled by the constant contributions coming 
from the terms proportional to c\ and C5. Up to 0(p 2 ), it is easy to work out from the Lagrangians 



(2) 

and C^Ljy the following expressions for the masses of protons, m p , and neutrons, m n : 



N 



m p = m N SBoCirh — ABoc^m , 

o 

m n = m N SBoCiih + ABqc^ui , (B.5) 

where rh = {m u + m<i)/2 is the average light quark mass and fh = {m u — m<i)/2 measures the strength of 
strong isospin breaking. When inserting the previous expressions in eq. flB.4| ) the difference I2 — Mjv 
will cancel those terms from c\ and C5 that are independent of the quantum pion fields or the external 
sources v, a, s and p. Indeed one can generalize this result for any operator, with n > 2, following 
similar arguments to the ones given in ref. . Among the building blocks for the construction of a chiral 
Lagrangian, only the scalar operator \+ gives rise to terms without pion or external source legs, 
leading to a constant matrix. Given x+i one can then build up operators O to be sandwiched between 
the t/j and ip fields which give rise to constant matrices without any field. In addition, one can include 
any number of covariant derivatives D^, since they count as 0(1) and the partial derivative d^, contained 
in D^, does not include any pion or external source in contrast to T^. However, we can integrate by 
parts in such a way, that all the <9 M act either to the left or to the right of the x+i as otherwise they 
will give rise to vertices with at least one pion or external leg. When all the Lorentz indices of these 
operators are contracted with the 7— matrices, the metric tensor and the totally antisymmetric tensor in 



(1 = 4 dimensions [34|, it follows from the application of the equations of motions and from the trivial 
result [d^jdu] = that the operators in the covariant derivatives just give rise to factors of m. These 
factors can be reabsorbed in the couplings of those operators that are only constructed from x+ ° r i n the 
vertices with one or more legs. As a result all the vertices without any legs are just constant and their 
finite parts are reabsorbed in the definition of the physical nucleon mass matrix Mat. The infinite parts 
that could arise because of possible scale dependent counterterms in front of such terms, needed for the 
renormalization, have to be considered jointly with pions loops, which obviously involve pion legs. 



C Application of the non-standard counting rules to in-medium pion 
propagation 

As an illustration we will discuss the application of the non-standard power counting in detail, as it 
appears for the in-medium pion-propagation. The leading contribution is C(p 3 ) and contains one gener- 
alized vertex, V p = 1 at lowest order, 5\ = 4, while V n = 0, n m = for m > 2 and Ib = I*b- Obviously 
p m = for all m since = 0. This case corresponds to figs.||a, where only one Fermi-sea line has been 
shown although it should be clear that one can substitute any thin line by a thick one including the 
appropriate factors as discussed in detail in appendix [A|, in sec. 4J and in ref. [|]|. The next-to-leading 
order is then 0(p 4 ) and there are various contributions: 

a) The same as before but with AIb = Ib — Ib = ^- This corresponds to fig.||b. Other diagrams 
similar to this one obtained just by exchanging the positions of the lines and the loop are not 
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shown. Instead of a loop one can also have an 0(p 3 ) irN counterterm without pion lines just for 
renormalization, as discussed at the end of sec. [b]. 

b) V p = 1, 5\ = 4, V,r = 0, ni = 1 and AIb = 0. This gives rise to the diagrams in fig.|5|c,d, where, 
as before, diagrams obtained under the exchange of the positions of lines and the loop or by the 
inclusion of any extra Fermi-sea insertion are not explicitly shown. The diagram of fig. ||d is zero 
because the Weinberg term with two pions involves the vertex [</>, <9 M </>] = 2i€ijk4>id fl (f)j, and hence 
one cannot have the same two flavors in one vertex. 

c) V p = 2, <5 1; 2 = 4, = 0, n m = m > 2 and AIb = 0. The diagrams are shown in figs.[|]e,f. As 
before diagrams obtained by exchanging lines and adding more Fermi-sea insertions are not shown. 
The diagram in fig. ^f vanishes because the new generalized vertex with only one pion is zero since 
there is a 75 matrix from the pion vertex and at most two additional gamma matrices. 

It is worthwhile to mention that in this example there is no contribution at next-to-leading order from 
the insertion of one operator in a generalized in-medium vertex, because then there are at least two 
pion lines involved, so that 5± = 5 and ri2 = 1. As a result such contributions only start at C(p 5 ). 

In order to calculate the diagram in fig.|5|a one has to make use of the full baryon propagator and not 
the heavy baryon (HB) expansion in eq. ( |2.1| ). We have calculated this diagram for the case of symmetric 
nuclear matter, where, as discussed above, the diagram fig.|||c vanishes as well. Then, up to 0(p 3 ) and 
for p = 0, we have instead of eqs. ( ^.5| ): 

—i 



Q 2 - Ml + 2m »&Q 2 (J( ^ Q) + x{ _^ Q)) 

|^ Uu)Q{k% - 2m N u)9{k F - M" / f " 4| Q J [ ( " ) +^ -«) = , (CI) 



where kj^(uj) is the maximum of k 2 F - 2m N u; and (Q 2 + 2m N u) 2 /4Q 2 . The function X(Q°, Q) is given 
by: 



J(Q°,Q) = V J" 



dp 



(2vr) 3 Q 2 + 2m N Q° - 2pQ 



k F {Q 2 + 2m N Q°) {Q 2 + 2m N Q ) 2 - Ak 2 F Q 2 (Q 2 + 2m N Q° - 2k F \Q\ 



,2 



16tt 2 Q 2 128tt 2 |Q| 3 6 (Q 2 + 2m N Q° + 2k F \Q\) 2 ' 

(C.2) 

where V indicates the principal value of the integral. The step functions in eq. ( |C.1| ) stem from the 
interplay of the second Fermi-sea insertion in fig. ||c with the pole from the free baryon propagator in 



fig.Hb as discussed in sec. 4.1 



D Wave function renormalization 



The pion fields are naturally normalized by imposing the canonical commutation relations at equal times 



y 



0. 



n(x) a ,IL(yf =i5(x-y)5 a P 



7T(xrMyf] = \u(x) a ,U(yf 



, 



(D.l) 
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where the indices a, f3 can refer to any of the tt°, ir + and tt~ fields and H(x) a is the conjugate momentum 
of the field ir(x) a . We now proceed to calculate the in-medium pion field and conjugate momentum from 
the quadratic term in the pion fields, f dxC^, of f dxC 

The in-medium free pion fields are given, at the classical level, by the equations of motion, 

S / dxZ^/84>k{y) = , 



already discussed in sec. These are fulfilled by writing the pion fields as: 

dQ 



7T + (x) 

ir~(x) 
ir°(x) 

where Q°(Q 2 ) = w+(Q 2 ), Q°(Q 2 ) 



(2tt) 3 
dQ 

(2vr) 3 
dP 

(2^)3 



a(Q)e- iQx + tf(Q)e iQx 
b{Q)e~^ x + a\Q)e iQx 
c(P)e- iPx +c t (P)e iP 



(D.2) 



-(Q 2 ) and P°(P 2 ) = wo(Q 2 ) are positive functions of the square 
of the three-momentum and correspond, respectively, to the 7r + , 7r~ and 7r° energies and are given by the 
corresponding spectral relation between energy and three- momentum, eqs. (|4,5|). Note that in eqs. ( p.2| ) 
we have space-time dependences of the form e~ l ® x and e l ® x . When applying eqs. ( |4.5[ ) to the latter case, 
we have to substitute Q° — > — Q° since these formulae were deduced taking a generic form e~ l ® x and 
allowing for both positive as well as negative values of Q° = uj. Taking this into account, together with 
the fact that the terms proportional to p in eqs. (|4.5| ) are odd in u, we can then combine the coefficients 
a(Q) and b(Q) and their conjugates in the 7r + and tt~ fields as shown in eqs. ( p.2| ). On the other hand, 
the in-medium contribution to the pion conjugate momenta ILi(x) are calculated by differentiating the 
quadratic piece J dxSC^ (given up to 0(p 5 ) and times a factor of i in eq. (4.3)) with respect to 



In order to deal with the presence of the non-localities induced by the baryon propagator Dq 1 (x,i/), we 
have worked in momentum space after calculating the functional derivative. Proceeding in this way, one 
finds for IT(x): 



Il(x)i = 4>(x)i 



{ 1 + f 2 ( C2 + C3 ~ } + ti e%k W 2 ^ 



9aP 



dQ Q 2 



2m N p J (2tt)3 2u, vac (Q 2 ) 2 



i(Q)e 



(D.3) 



where we have denoted the vacuum energy by w vac (Q 2 ) = y M 2 + Q 2 . One can easily calculate the con- 
jugate momenta of the fields 7r + (x), tt~(x) and tt°(x) from eq. ( p.3j ). Taking also into account eqs. ( p.2j ) 
one has: 



U(x)- 



u( x y 



p_ 

2p' 



dQ 

(2tt)' 



dQ 

(27T)* 



• P +( \ 

+ l 2p % W 



LO-{Q 2 )b{Q)e~ lQx - uj + {Q 2 )a\QyQ x 



1 + 



4 P , 

J2 ( c 2 + C3 



(Q 2 )a(Q)e- iQx - uj _(Q 2 )&t(Q)e iQa 



A m 2 

8m N cJ vac (Q 2 ) 2 

+ /2 ^ 2+C3 g miva;vac(Q2)2 
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n(x)° 



dP 



■M? 2 ) c(P)e 



-iPx 



c f (P)e 



iPx 



1 , ^ I 

1 + -p I C2 + C 3 



4 



M 2 



8m at w vac (P 2 ) 2 



(D.4) 



We now proceed to the quantization of the pion fields eqs. ( p.2[ ) by imposing the canonical commutation 
relations given in eqs. (DA) for equal times. We denote by d a (Q) a generic coefficient in eqs. (p.2|) and 
by <ia(Q) its complex conjugate, with a = referring to a 7r° and a = +(— ) to a 7r + (7r~) state. We 
quantize them by imposing the following commutations relations: 



[d a (Q),dp(P)} = (27rfN a (Q 2 )5(Q-P)6 af , , 



[d Q (Q),^(P)]= [4(Q),4(P) 



. 



We now prove that it is possible to satisfy eqs. ( p.l| ) by an appropriate choice of iV a (Q 2 ). 
iVo(Q 2 ) is fixed by the conditions x° = y° and 7r°(a:), IT (y) = i<5(x — y) with the result: 



(D.5) 



N (Q 2 



2u; ^(Q 2 ) 



where 



Z(Q 2 ) = 1 + -£ ( c 2 + c 3 



f 2 



Ml 



(D.6) 



(D.7) 



8mAT Wvac(Q 2 ) 2 7 ' 

The other commutation relations involving the 7r°(x) fields are then trivially satisfied. In the case 
of the charged pion fields, the commutation relations at equal times read: [tt + (x), tt~ (y)] = and 

7r + (x),n + (y) = ir~(x), H~(y) = iS(x — y). These conditions imply the following equations: 



iV + (Q 2 )=Ar_(Q 2 ) , 

o;+(Q 2 )iV + (Q 2 )+a;_(Q 2 )iV_(Q 2 ) = Z(Q 2 )- 1 



Thus: 



iV + (Q 2 ) = iV_(Q 2 ) 



1 



(D. 



(D.9) 



K(Q 2 )+ W _(Q 2 )]Z(Q 2 ) • 

We are now in the position to discuss the wave function renormalization. We will normalize our pion 
states with definite three-momentum p, |7r Q (p)), as: 



(vr Q (q)|vr Q (p)) = 2 WQ (q 2 )(2vr) 3 J(p-q) 



(D.10) 



Here we have chosen the functions in front of 5(p — q) as in the vacuum, so that the calculated matrix 
elements do not change with respect to the vacuum ones in the limit kp — > 0. On the other hand, we 
define the one-pion state in terms of the creation operators as: 



K(p)> = e a ( P 2 )4(p)l^> . 



(D.ll) 



Imposing eq. ( p,10[ ), together with eqs. ( |D.5| ), one has (except for a multiplicative phase that can always 
be chosen to be one): 



e Q ( P 2 ) 



'2a^(p2) 
N a (p 2 ) • 



(D.12) 
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In the calculation of a matrix element of a certain process with asymptotic pions, the latter are annihilated 
or created by the corresponding pion fields. In this way, taking into account eqs. flD,2D , ( D.5 ), ( D.ll ) and 

flrnD, one finds that: 



(n\ir a (x)\7r p ) = ^/2u a N a {p 2 )S aP e 



ipx 



(D.13) 



with Z a (p 2 ) the wave function renormalization of the pion states of type a in the nuclear medium. As a 

result we can proceed directly by annihilating and creating asymptotic pion fields from the fields ir a (x), 

_ i 

neglecting the factors Z a 2 during the calculation, but multiplying in the end by as many of such factors 
as there are asymptotic pion states, 

l[Z a (p 2 )-l (D.14) 



The explicit expressions for Z a (p 2 ) are: 



z + (p 2 ) 

Z-(P 2 ) 
Z,{p 2 ) 



1+^ 

p 
p 



1 + 



p 



C2 + C 3 
C2 + C 3 
C2 + C 3 



9 2 A 



M, 



8m N w vac (p 2 ) 2 
Ml 



9aPP 2 



2pu- 



2, -VacP 2 ) 



2/ 2 ^vac(p 2 ) 3 



9\ 



8m N w vac (p 2 ) 2 

9\ Ml 
8m N u; vac (p 2 ) 2 



+ 



9aPP 2 



2/ 2 o, vac (p 2 ) 2/ 2 Wvac (p 2 )3 



(D.15) 



It is worthwhile to point out that the wave function renormalization in eqs. ( D.15| ) cannot be obtained 
from eqs. ( ^4.5| ) just by calculating their derivatives with respect to uj 2 at a value determined by the 
equations of motion (for a given three momentum). The reason is the non-local character of 5C^, see 
eq. ( |4.3| ), This implies the presence of extra factors of energy, to, that stem from the propagators of the 
nucleons and have nothing to do with time derivatives of the pion fields, as in standard local quantum 
field theory. 

We do not discuss the vacuum contribution to the wave function renormalization, since we take directly 
the final 0(p A ) results from || and since the effects of mixing the vacuum and in- medium wave function 
renormalizations are of higher orders, (D(p®) in the standard counting. In the latter reference the vacuum 
Green functions of the external sources were directly calculated in terms of the generating functional up 
to one pion loop. One has to stress that the calculation of the Green functions of the external sources 
connects directly with real experiments where the nuclear medium is tested by means of external probes, 
e.g. muon capture, neutrino scattering, etc. This is so because the vector source v can be related to a 
photon and the axial-vector a to the gauge bosons of the weak interactions. In this way one can 
incorporate electromagnetic and semileptonic weak interactions in the formalism, see e.g. ref. p0|. 
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